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1. Contents of the lectures 

The lectures are organized as follows, in Sections [2] and [3] the properties 
of QCD matter and the use of heavy-ions to study these extreme conditions 
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are briefly reviewed; Section [J] presents the general formalism to describe 
highly energetic particles propagating in matter which is then applied to the 
case of medium-induced gluon radiation in Section [5j Section [6] summarizes 
how this formalism is implemented in actual phenomenological calculations 
and gives a brief review of the comparison with experimental data; Section 
[7] presents some new developments on the calculation of the transport co- 
efficient using the Anti-de-Sitter/Conformal-Field-Theory correspondence; 
finally, some conclusions are presented. Readers interested on the formalism 
of medium-induced gluon radiation could directly jump to Section |U The 
derivation of the medium-induced gluon radiation presented here is new and 
entirely based on resummation of the multiple scattering diagrams which 
we find more intuitive. 



2. Heavy-ion collisions and extreme QCD matter 

QCD is the theory which describe the strong interaction. Its asymp- 
totic states are hadrons which are colorless objects composed by quarks and 
gluons, the degrees of freedom of the QCD lagrangian. Confinement is a 
property of the strong interaction which forbids the existence of asymptotic 
colored states in normal conditions of temperature and densities. Quarks 
are fermions with fractional electric charge and very different masses. Light 
quark masses are O(10 MeV) for up and down and 0(100 MeV) for strange. 
Heavy quark masses are ~1.5 GeV for charm, ~5 GeV for beauty and ~175 
GeV for top. The masses of the hadrons formed by light quarks is much 
larger than the sum of the individual masses of their constituents, indicating 
a large dynamical origin for the first. The spectroscopy of both light and 
heavy mesons and baryons is an interesting area to study the QCD dynam- 
ics. A main observation is that the chiral symmetry - an exact symmetry of 
the lagrangian when the quark masses are neglected - is not observed. This 
symmetry breaking is not realized by the presence of a new particle as the 
Higgs boson in the EW sector, but by the structure of the QCD vacuum, in 
which the symmetry is spontaneously broken. In this case, the associated 
Goldstone bosons are identified to be the pions and kaons, which are almost 
massless when their mass is compared with that of other hadrons. 

So, two basic properties of the QCD vacuum are the confinement of 
quarks and gluons and the chiral symmetry breaking. One relevant question 
is, then, is there a regime where these symmetries are restored? 

Another essential property of QCD is asymptotic freedom, the fact that 
the interaction disappears at small distances/large scales. As a consequence, 
QCD medium at asymptotic temperatures is predicted to be a gas of free 
quarks and gluons with a restoration of the symmetries mentioned above. 
For a null net baryon number in the medium, the transition temperature 
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Fig. 1. Left: Entropy, s = e + p in units of T 4 versus temperature computed in 
lattice. Right: Trace anomaly in units of T 4 as a function of the temperature. 
Figures from [3] 



is found to be in the range T c ~ 175 - 190 MeV - see e.g. [JJ for a recent 
discussion. Most of the theoretical information we have about the transition 
region comes from lattice QCD calculations where different quantities can 
be studied^ such as the equation of state, the chiral condensate or the free 
energy. 

A first example of this collective behavior is the equation of state. The 
pressure or the energy density measured in units of T 4 of an ideal gas are, 
according to the Stefan-Boltman law, proportional to the number of degrees 
of freedom in the system: so, for a free gas of pions this quantity is NT* = 3, 
while for a free gas of quarks and gluons this quantity is much larger, JV? - = 
2x2x3, N% * = 2x8, counting spin, color and (two) flavor states. This 
different behavior is observed in lattice calculations where a jump at the 
transition temperature, T c , appears - see Fig. [TJ Interestingly, the lattice 
results also signal to a significant departure of the ideal gas behavior, e = 2>p 
at temperatures close to T c - see Fig. [TJ 

In order to understand the nature of the transition order parameters 
are needed. In QCD with massless quarks, the chiral condensate is the 
order parameter, while in the infinite mass limit the order parameter is the 
Polyakov Loop. The order of the transition depend on the actual value of 
the quark masses and most calculations agree in the absence of a real phase 
transition at zero baryochemical potential - the transition would just be a 
rapid cross-over. 

The behavior of the potential between heavy quarks are also of interest 
for phenomenological applications. Although some discussion about the 
precise meaning and the definition of a potential exists in the case of a hot 
medium a screening leading to a non-confining potential is expected to 
appear at some point above the phase transition and, correspondingly, the 



See e.g. [2] for a description of the QCD thermodynamics as studied by lattice. 
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heavy quark bound states would cease to exist if this temperature is reached. 
A prediction of the presence of a deconfined medium in heavy ion collisions 
is then the disappearance of heavy-quark bound states (quarkonium) and 
in particular the J/\l/. 



The description of hadronic collisions in terms of the formation of a 
thermal state, which then evolves according to a hydrodynamical behavior, 
has an old story. In the 50's, Landau [5] proposed a model in which a 
transient thermal state - little fireball - is formed by energy deposition in a 
small slice of space, of the size of the Lorentz-contracted nuclei, at the center 
of mass of the collision. The subsequent hydrodynamical evolution produce 
an expansion and cooling of the system up to a freeze-out temperature 
T ~ when the formed hadrons free stream to the detectors. A rapidity- 
independent version of this model was proposed by Bjorken [6] - see e.g. [7] 
for a review on hydrodynamical models in heavy-ion collisions. Although 
more sophisticated implementations are used nowadays, very demanding 
from the computational point of view, for the hydrodynamical description 
of nuclear collisions, the Bjorken model estimate gives a convenient and 
simple parametrization of some medium properties, as a function of the 
proper time r, as the energy density or the temperature of the created 
medium 



Present facilities to study this active field of research - see [8] for a recent 
review on the experimental situation - are the CERN SPS, with fixed target 
collisions of different nuclei up to ~ 20 GeV/A; the RHIC at Brookhaven, 
a dedicated collider which reaches much larger energies, y/s =200 GeV/A; 
and in the near future the LHC at CERN, with an important program 
of nuclear collisions which will provide a spectacular jump in the collision 
energy upto \fs = 5500 GeV/A. Here, the quoted figures refer to the total 
available energy in center of mass of the collision divided by the number 
of nucleons in the nuclei - energy per nucleon - so that it can be directly 
compared to a corresponding proton-proton collision energy. 

The lifetime of the hot medium produced in heavy ion collisions is so 
small, of the order of its own transverse size, that only indirect probes are 
available to characterize its properties. Here we will only focus on hard 
probes and specially on the study of the high transverse momentum part 
of the spectrum of produced particles. Jet quenching is the generic name 
under which the corresponding effects on high-pj_ particles are known. The 
original proposal [9j, a suppression of the the large transverse momentum 
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yields due to the energy lost by interaction of the fast partons with the 
medium has been observed at RHIC [10] and partially also at the SPS [TT] . 



3. Hard probes to heavy-ion collisions 

Hard processes are those involving large momentum exchanges. Asymp- 
totic freedom allows, in these conditions, to perform calculations by a per- 
turbative expansion in terms of a s (Q 2 ), where Q is the large scale of the 
process, as, e.g. a large mass or a large transverse momentum of the pro- 
duced particles. This perturbative expansion is, however, computed in terms 
of quarks and gluons, the degrees of freedom of the QCD lagrangian, and 
some extra information is usually needed to connect it with the initial and 
final hadrons appearing as asymptotic states. A crucial simplification of 
the problem is possible thanks to the factorization theorems of QCD which 
allow to separate long- and short-distance contributions to the cross section 
in the forrrQ 

a AB-+h = /a(xi> Q 2 ) 3 Q 2 ) 3 a ( XliX2j Q2 ) g, D ^ h {z, Q 2 ) , (2) 

where the short-distance perturbative cross section, a(xi, X2, Q 2 ), is com- 
putable in powers of a s (Q 2 ) and the long-distance terms are non-perturbative 
quantities involving scales O(Aqcd)- More specifically, the proton/nuclear 
parton distribution functions (PDF), /^(x,Q 2 ), encode the partonic struc- 
ture of the colliding objects and can be interpreted as the probability to find 
a parton (quark or gluon) inside the hadron with fraction of momentum x; 
the fragmentation functions (FF), D(z,Q 2 ), describe the hadronization of 
the parton i into a final hadron h with a fraction of momentum z. 

Although the PDFs and the FFs are non-perturbative quantities, their 
evolution in Q 2 can be computed by the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) equations [12]. This equations resum higher orders in a s 
which, due to enhanced contributions in some regions of phase-space, can- 
not be neglected. For example, we will see Section [5] that the spectrum of 
gluons radiated with energy fraction z out of a quark produced in a hard 
process is 

dl a s 1 a s l Cr 



dzdt 2vr 2 t P(z) ~ 7r 2 t z (3) 

with t the invariant mass of the produced pair. This radiation is divergent 
in the infrared and collinear limits giving large 0(a s \ogt) terms to the cross 
section which need to be taken into account. These terms are resummed by 



2 For simplicity we have made all the scales the same, although, strictly speaking, the 
factorization, renormalization and fragmentation scales can be different. 
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the DGLAP evolution equation^] 



dD(x,t) 
dlogt 



L 



ld ^P{z)D( X -A. 
z 2vr K 1 \z J 



(4) 



where D(x,t) can be the PDFs or the FFs and P(z) are the splitting func- 
tions. At the lowest order (LO) these equations admit a probabilistic in- 
terpretation in which a parton shower is formed by subsequent branching 
of the initial quark or gluon. We will come back to this interpretation in 
Section El 

These equations provide an excellent description of experimental data 
as deep inelastic lepton-proton scattering (DIS) or jets. 



The hard process described by the perturbative cross section a(x%, x%, Q 2 ) 
in ([2|) takes place in a very short time, 0(1/Q), and is basically unchanged 
in nuclear collisions. The interest of hard processes as probes of the medium 
is that the long distance parts are modified when they communicate with 
the extension of the medium. Measuring these modifications allow for a 
characterization of the medium properties - see e.g. [13] for a recent re- 
view. 

A conceptually simple example is the J/^, whose production cross sec- 
tion can be written as 



a hh^j/* = fi{xuQ 2 ) g, /.( X2 ,Q2) ® a i ^(x 1 ,x 2 ,Q 2 )(0([cc] -> J/*)) , 



where now (C([cc] — * J/^f)) describes the hadronization of a cc pair in a 
given state (for example a color octet) into a final J/^f. In the case that the 
pair is produced inside a hot medium this long-distance part is modified. 
Actually, the potential between the pair is screened in a hot medium and 
the hadron is dissolved, making (0([cc] — ► J/^f)) — > 0. The experimental 
observation of this effect is a disappearance of the J/^f in nuclear collisions 
[15] . This suppression has been discovered in experiments at the CERN SPS 
|16j and then also at RHIC [TTj. The interpretation of the experimental data 
is, however, not simple; one of the reasons being the lack of good theoretical 
control over the modification of a purely non-perturbative quantity as the 
hadronization of a cc pair into a charmonium state. 

From the computational point of view, a theoretically simpler case is 
the modification of the evolution of fragmentation functions of high-pr par- 
ticles due to the presence of a dense or finite-temperature medium. Here, 

3 To simplify the notation, we include only one flavor, see, e.g., [13] for a complete 
description. 



3.1. Hard processes as probes of the medium 
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highly energetic partons, produced in a hard process, propagate through the 
produced matter, loosing energy by medium-induced gluon radiation. This 
phenomenon is generically known as jet quenching and its implementation 
into the perturbative formula (|2|) will be discussed in the next sections. 

3.2. Nuclear parton distribution functions 

A good knowledge of the PDFs is essential in any calculation of hard 
processes. The usual way of obtaining these distributions is by a global fit 
of data on different hard processes (mainly deep inelastic scattering, DIS) to 
obtain a set of parameters for the initial, non-perturbative, input f(x, Qq) 
to be evolved by DGLAP equations |12j . 

In the nuclear case, the initial condition, Ja(x, Qq), is modified compared 
to the proton. Moreover, at small enough x, non-linear corrections to the 
evolution equations [18j QjJJ [201 [21] are expected to become relevant. Global 
DGLAP analyses, paralleling those for free protons are available |22 } I23| , l24" t 
[221 [23 123 ES]- These studies fit the available data on DIS and Drell-Yan 
with nuclei providing the needed benchmark for additional mechanisms. 

The DGLAP analyses of nuclear PDFs from [27] is shown in Fig. [21 
including the corresponding error estimates. An important issue, partially 
visible in Fig. [21 is that present nuclear DIS and DY data can only constrain 
the distributions for x > 0.01 in the perturbative region. By chance, this 
region covers most of the RHIC kinematics, so that, the description of e.g. 
J/^-suppression or inclusive particle production in dAu collisions as given 
by the nuclear PDFs can be taken as a check of universality of these distri- 
butions. These checks present a quite reasonable agreement with data |29| . 
but some extra suppression for the inclusive yields at forward rapidities is 
probably present. The strong gluon shadowing plotted in Fig. [2] improves 
the situation at forward rapidities without worsening the fit of DIS or DY 
data - x 2 /dof < 1. Whether a DGLAP analysis can accommodate all sets 
of data is an open question, but the finding in Ref. [2Zj are encouraging. A 
suppression at forward rapidities was also predicted in terms of saturation 
of partonic densities [30] . 

4. Particle propagation in matter 

To describe the jet quenching phenomenon, we start from a high-pj_ 
quark or gluon produced in an elementary hard collision which subsequently 
interacts with the surrounding matter. So, the first question we need to ad- 
dress is how a highly energetic particle propagates through a dense medium. 
In this section we will present the most widely used formalism based on a 
semiclassical approach in which the changes in the medium configuration 
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A = 208, Ql = 1.69 GeV 2 
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Fig. 2. Ratios of nuclear to free proton PDFs for different flavors at the initial scale 
Qo=l-69 GeV 2 from [STJ with error estimates. The green line in the gluon panel is 
an attempt to check the strongest gluon shadowing supported by present data. 



due to the passage of the fast particle - recoil - are neglected. In this 
approximation the medium can be considered as a background field. 
A convenient formulation of the problem is in terms of Wilson lines 



W(x) =Vexp 



ig / dx+A-(x+, x) 



(6) 



describing the propagation of a particle through a medium field A_(x + x). 
Its origin will be explained in the next subsection. Here we introduce the 
light cone variables 

1 , . 1 . 

X± = — (SCO ± X3) P± = —7= (P0 ± P3) 



So that the scalar product is 

p ■ x = p + x- + p^x + — p± • Xj_ 

and the rapidity 

V 



2 



P0 +P3 



P0 ~P3 



2 



P- 



(7) 

(8) 
(9) 



4-1- Wilson lines, eikonal approximation 

A simple derivation of the Wilson line is obtained in terms of multiple 
scatterings, providing a clear physical picture of the eikonal propagation. 
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X\ x 2 x 3 x 4 x n 
Fig. 3. A multiple scattering eikonal trajectory 



Consider the diagram of Fig. [31 where static centers of scattering are placed 
at x\, X2, ... x n . Let us fix that the quark is moving in the positive X3 
direction, i.e. the large component of the momentum is p+. 
The contribution to the S-matrix of one scattering is 



Si(p',p) = / d A x S p '~ p) - X u(p') igA a Jx)T a ^ u(p) 



(10) 



Taking the eikonal limit, p ~ p', \'}ZxU X (p)^u x {p) = 2p fl and p^Af^ ~ 
2p + A°_ . In order to proceed, we will assume that the fields have a small 
dependence on the small coordinate x_: due to the Lorenz contraction, the 
medium can be seen as a small sheet in this coordinate. Putting all together 
one obtains 

S 1 (p',p)^2ir5{p' + -p + )2p + J dkj_e~ <Xj - (p J-~ p - L ) ig J dx + A_(x + ,x ± 

(11) 

where we have singled out with brackets the contribution of the field which 
will exponentiate to give the Wilson line and the color matrix has been 
omitted for clarity. 

The contribution with two scatterings is given by 



S2(p',p) 



^ I d A x 1 d A x 2 e^i-^V^'-^^W) igA a u \{ Xl )T a ^^ 

(27T) 4 V" 1 



X I 



^-igA^(x 2 )T a ^ u{p) 



pj + ie"^^ 2 



(12) 



4 Note that we neglect the p- component. In the eikonal approximation p'_ << |pj_|, 
however, the phase factor ip_x+ is potentially enhanced by the medium lenght. Thus, 
we are implicitilly assuming that the medium is smaller than the coherence length 
P-L ~ fi 2 L/2p+ << 1. In the next section we will see how this assumption can be 
relaxed. 
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Applying the Dirac equation one can write in the eikonal limit u(p)jui i>il^2 u {p') 
(2p+) 2 5mi-^ 2 - giving 



e i(pi—p)-xi+i(p' —pi)-xn 

ct i xid' i X2 A-(xi)A-(x 2 ), 



Pi + It 



( 13 ) 

where again the color matrices have been omitted as will be omitted in the 
following. In the high-energy limit, the integrals in d 4 pi can be performed: 



e i(asi+-as2+)pi- 2iri 

dpi-- — — = -Q(x 2 + - xi + )- — 

2pi+Pi- + k 2p + 

dp 1+ e iixi — X2 - )pi + = 2ir5(xi- - X2-) 
J cPp 1±e - i( - x ^- x ^P^ = (2vr) 2 5(x 1± -x 2 ±) 
Doing the ^--integration as in the case of one scattering one obtains 



(14) 
(15) 
(16) 



1 2 



S 2 (p',p) ~2x5(p' + -p + )2p + J dx ± e- ix ^- p±) ^V ig J dx + A^{x + ,x ± 

(17) 

where again the second term in the expansion of ([6]) appears. In order to 
obtain (fT71) we have used that 



/dxidx2-.-dx n Q(x2—xi)....@(x n —x n -i)A(xi)...A(x n ) = —V f dxA(x) 
n! [J 

(18) 

where V means path-ordering of the fields A(x). 

The generalization to n-scattering centers follows the same lines and the 
corresponding contribution is given by 

S n (p',p) ~ 2ir5(p' + - P+ )2p + J d x± e-^(P'x-P^lp ig J dx + A„(x + ,x ± ) 

(19) 

So that the total S-matrix is 

oo . 
S(P,P) = Y, S n(P^P) - 2 ^(P+ -P+)2P+ / dx ± e- ix ^-P^(x ± ), 



n=0 



with W(xj_) given by eq. Q. Similar expressions are obtained in the case 
of a fast gluon of momentum k with the only replacement of the external 
fields in ([6]), now in the adjoint representation - see e.g. [31]. In order to see 
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this, the external fields coupled to three gluon vertices (four gluon vertex 
are subleading in energy) are written in the adjoint representation 

A A = A a (T%) bc = -iA a f abc (21) 

The sum of gluon polarization vectors can be written as 

2 

m^k u k fl m u 



ww m ■ k m ■ k 



i=i 



with m a light-like vector with non-zero minus component, m = (0,1,0,0) 
and the polarization vectors defined by e • k = e • m = 0; e 2 = — 1 

HV = (°>^-^>t(i)±) (23) 



and = (1)0), £(2)i = (0,1). This prescription allows to simplify the 

structure of the gluon line since, due to gauge invariance, only the first term 
in eq. (122f) contributes in the high-energy limit when the propagator is in 
between two three-gluon vertices. So that, one can write for each gluon 
propagator 

i=i 

and now the contribution from each vertex Vp, va can be simplified as 

e1 if) (k')V^ a (-k',k,k' - k)A^e u {i) (k) ~ -i2gk + A A „5u> (25) 



which has a structure similar to the quark case allowing for a resummation 
to obtain a Wilson line now with fields in the adjoint representation. We 
will denote this gluon line with a superscript, W a (~k±). 

Eq. (|20|) describes the scattering of a quark (or gluon, with the changes 
explained above) in a medium with a given field configuration, determined 
by the scattering centers. The effect of the medium is to induce color rota- 
tion at each scattering center without changing the helicity (polarization) 
of the quark (gluon), which loses a negligible amount of energy and prop- 
agates in straight lines. In order to have the final answer for a physical 
cross section, these medium configurations need to be averaged within an 
ensemble describing the medium. This issue will be discussed later. 
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4-2. Relaxing the eikonal approximation 

In some cases, the restrictions applied in the above formulation need to 
be relaxed to allow small changes in the transverse position of the prop- 
agating particle. This is the case, for instance, of the medium-induced 
gluon radiation, where the gluon position follows Brownian motion in the 
transverse plane. The eikonal Wilson line is now replaced by the propagator 



G(b,a) = J Dr(x+) exp < i^- j dx- 



dr 

dx + 



W(r) 



(26) 



This change can be derived in the multiple soft scattering presented in the 
previous section. In order to do that, we have to keep the subleading p 2 ^ 
terms in the poles of the propagators. So, now the integration in pi- reads 



dp- 



Jp-(x i+ -x {i+1)+ ) 



2vr 



e(x 



(t+l)+ ~~ x i+) e 



2p- 



;((«»+ -*(»+!)+) 



2p+p_ - p 2 ± + ie 2p + 
and instead of (fl6j) the integration in pi± is gaussian, giving 



(27) 



(27T) 2 
P± 

2m(x i+ - x (i+ i) + ) 



»+ :r (i+l)+) e -«Pi±(x i± -x (i+1)± ) 



cxp 



.P+ (Xi± - X( i+ i)j 



2 Xi 



(28) 



Eq. (|28p is the Feymann propagator of a free particle that propagates in the 
transverse plane from Xij_ at time Xi + to X(j +1 )j_ at time (see e.g. 

[33] ) Go(x(i +1 )i_ ~~ x i±! x {i+i)+ — x i+)- Thus, Eq. ([28]) can be expressed as 



G , o(x( i+ i)_ L - x i± ; - x i+ ) = j T>xj_(x p ) exp j i^- J dx A 



dx-x 
dx + 



(29) 



where the paths x_|_(x+) connect the two endpoints of the propagator. 
The scattering matrix now reads 

oo „ n 

S(p',p) = 2tt5(p' + -p + )2p + ^2 / 'PY[dx i+ dxi±igA(xi + ,Xi±) x 

n=0 J i=0 

xGo(Ax (i+1)1 ; Ax( i+1)+ )ic/A(x (i+1 ) + ,x (i+1) _ L ) (30) 

where Ax( i+1)1 = X( i+1 )j_ — xll, Ax (l+1)+ = x^ +1 - )+ -x i+ . This expression 
can be reorganized as 



S(p',p) = 2tt5(p' + -p + )2p + I dx ± e 



^-ixxCp', -pj 
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with W(x±) given by ©. 

Eq. ([3T]) describes the propagation of a highly energetic particle through 
a medium when changes in the transverse position are allowed. The prop- 
agation takes phases both from the motion in the transverse plane and the 
color rotation by interaction with the medium field. 



4-3. The medium averages 

The S-matrices derived in the previous sections are valid for a given 
configuration of the fields and should be averaged over the proper ensemble 
of the medium field configurations. Several prescriptions have been used in 
the literature. Here we will present two of them used for the calculation 
of the jet quenching. Notice first that any physical quantity must contain, 
at least, the medium average of two Wilson lines since the average is done 
at the level of the cross section where only colorless states are allowed. In 
general, we will be interested in quantities like 



^Tr(W\xL ± )W(y ± )). 



(32) 



Where the trace and the 1/N term correspond to average the initial color 
indices - correspondingly, a factor 1 / (N 2 — 1) would appear in case of Wilson 
lines in the adjoint representation. 

Two main approximations to the averages (I32p are used in jet quenching 
phenomenology, the multiple soft scattering approximation and the opacity 
expansion. Both can be understood in the multiple scattering picture we 
are employing. The main assumption is that the centers of scattering are 
independent, i.e. no color flow appears between scattering centers separated 
more than a distance A ~ l//i, \x being a typical scale in the medium as the 
Debye screening length. So, we want to calculate 

^Tr(I^t(x ± )^(y ± )> = l T r^exp{-^ J dx + A ] _ (x+, x ± )} x 

x exp{ig / dx+A^(x + ,y_ L )}) (33) 



Expanding the exponents and taking the contribution from one scattering 
center, the leading contribution is quadratic in the fields - the linear con- 
tributions cancels due to the color trace 



v(l + \{ig) 2 



dx + A Jf _(x + ,x±) 



dx+A-(x+,y±) 
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(34) 



These terms have a clear diagramatic interpretation, they correspond to 
the two-scattering diagrams represented in Fig. 21 The first two do not 
resolve the transverse dimension of the dipole and are sometimes called 
contact terms. Writing the field as the Fourier transform of a field amplitude 
produced by a scattering center at position (a;j + ,X|j_) 

gT b A b _(x + , X± ) = J ^ e ^-^T b a b _( q )5(x + ,x i+ ), (35) 

the medium averages (...) are done by integrating in the longitudinal Xi + 
and transverse coordinates x n j_ of the scattering centers and adding color 
configurations. For example, the third term in (|34p gives 

j dx + dx, l+ dx i± J ^2 (2^)2 e-^- x '^V^- x '^ x 
(qi)a-(q 2 )5(x+ - x l+ ) = J dx+ J ^_ e ^~ x ^ |a(q)| 2 (36) 

Putting all together, eq. (|34p can be written in terms of the cross section 
of a dipole, where the quark and the antiquark are located at y± and x_|_ 
respectively, with one of the scattering centers of the medium. 

1 " °F [ 4% Kq)! 2 (1 " e^- x ^ q ) = 1 - la(y ± - xj.) (37) 



xa_ 



(2tt) 2 1 v ' 2 

The factor Cp comes from the color average and the traces over the matrices 

1 1\t2 _ 1 A 

— YTiT a T b = C F = — = - (38) 

2N 3 y ' 

a,b 

The corresponding factor for a gluon would be Ca = N = 3. 

There is a subtlety in the derivation we have just presented. Strictly 
speaking, the contact terms are not included in the multiple scattering 
derivation of the Wilson lines made before, where each scattering center 
contributes with a factor A- and path ordering appears. What we have 
implicitly assumed is that scattering centers separated by distances smaller 
than A ~ \j \x resolve color correlations while those separated by more than 
A do not. This is usually implemented by considering the contribution of 
one scattering center by defining a density n(x + ) = — x, i+) an< ^ 

taking the continuum by an integration in x+ 

^Tr(^t(x i .)VF(y ± )) 1 scatt ~ 1 - ^ J dx + n{x + )a{y ± - x ± ) (39) 
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Fig. 4. Different contributions to the dipole cross section 



-^Tr(WHx ± )W(y ± )) ~ exp ^~ J dx + n(x+)a(y ± ~ x±)j 



Eq. (|39p is the first order in an opacity expansion of the medium, the sum 
of all orders exponentiate and the average (|33p can be written as 

(40) 

Eqs. (|39p and (|40p are the two main medium averages used in the literature 
of jet quenching. In order to proceed, the functional form of the dipole 
cross section needs to be specified. In the opacity expansion a Yukawa-type 
elastic scattering center with Debye screening mass \jl is usually taken in 

(EZD 

|a(q)l 2 = - ( 2V (4i) 

When the number of scattering centers is very large, all of them need to be 
resummed and the first orders of the opacity are not enough. In this condi- 
tions, it is convenient to take the dipole cross section at leading logarithmic 
accuracy [M] and write the small distance component of the cross section 

a(r) ~ Cr 2 (42) 

The proportionality factor C with the squared dipole size is usually taken to 
be constant and defines the transport coefficient q(£) = 2y/2n(£)C, encoding 
all the information about the dynamical properties of the medium. This is 
the main parameter to be determined by fits to experimental data and to 
be compared with_theoretical calculations. The Wilson line averages define 
this parameter 



N 2 



^Tr(W A t(xi ± )W A (y ± )) ~ exp j-^-L J dx + q(x + )(x ± - y ± ) 2 | . 

(43) 



5 The factor y/2 is inciuded here as the transport coefficient is usuaify defined in ordi- 
nary coordenates, where the iongitudinai distance for a £ = z ~ x+/^/2 
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Here we have considered a gluon instead of a quark as the fast particle 
traversing the medium. This is the most relevant quantity as we will see in 
the next sections. 

It is worth noticing that similar medium averages are used in the study 
of deep inelastic scattering (DIS) with nuclei, where a phenomenon of satu- 
ration of partonic densities is expected. In this case, the convention is to call 
the relevant parameter saturation scale Q sa ,t and the Wilson line average is 
given by 

^Tr(Wt( x± ) W(y±) ) ~ exp |-Ig2 at(x± _ y±) 2 j (44) 

Now the saturation scale contains not only the dynamical information about 
the target but also the geometry. A simple model for this saturation scale 
gives [35] 

<2sat = N2 _ 1 R 2 xG ( X ' V(X± " y±) 2 ) (45) 

Where the factor A 1 ^ comes from the proportionality with the nuclear 
medium - the radius for a nucleus is approximately Ra — -Ro^ 1 '' 3 ) Rq — 1-2 
fm. In (|45j> . xG{x, Q 2 ) is the gluon distribution of a nucleon inside a nucleus 
giving an energy dependence to the saturation scale which, as we have said, 
is usually neglected for the case of q. 



4-4- A simple application: the dipole model 

The above rules allow us to write the cross section for the scattering of 
a (colorless) dipole of size r 2 = x 2 
scattering) amplitude as iT = 1 - 
a = 2ImT 



with a medium. Taking the (forward 
S and applying the optical theorem, 



a 



-Tr{WHx ± )W(0)) 



1 



exp 



.in 2 

a Vsat 



(46) 



This dipole cross section finds a direct application in the calculation 
of the deep inelastic scattering at small values of the Bjorken variable x. 
The dipole model tells us that in these conditions, the collision can be 
seen as a two-step process in which the virtual photon from the lepton 
splits into a qq dipole and then this fluctuation scatters with the hadronic 
target. The problem now reduces to compute the photon-dipole vertex ^j? L 
for a transversely or longitudinally polarized photon (see e.g. [36]) which 
factorizes to give the cross section for the virtual photon-hadron as 



a, 



Z L (x,Q ) 



d 2 r 



dz \ty 



7* 1 2 dip 



T,L\ 



o- + . 



(47) 
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Eq. (|47jl is extensively used in the phenomenology of DIS to include satura- 
tion effects - simpler to include in configuration space as we do here. As the 
photon wave function is known, the problem reduces, in general, to com- 
pute the dipole cross section for different medium averages and including 
different resummation techniques for the non- linear terms. 



4-5. Relation with the momentum broadening 

From the parton S matrix one can readily compute the momentum 
broadening. After passing though the medium, the particle distribution 
is given by 



-^-f- oc J dp' + 5(p' 2 -m 2 )^Tr(jS(p',p)[ 



d 2 p' ± 

And the transverse momentum broadening is given by 
/ 2\ 1 f ,2 / /2 dN 

From Eq. pip we find 

AA(pi)oc j dp' ± J dx ± dx' ± p' ± 2 e ip '^-^ 



(48) 



(49) 



VrVr' 



x exp < i- 



dx^ 



dr 

dx-i- 



dr' 

dx-L 



— Tr<VVW r )) (50) 



A simple manipulation leads to 
x^Ti(w\r')W{r) 



dx.1 



dr 

dx-L- 



dr' 

dx+ 



(51) 



In the high energy limit, p + » fj,, we can approximate the path integral by 
the classical path with smallest action, dr/dx+ = 0. Performing a similar 
manipulation for the normalization we obtain 



(pi) 



Ti(Wt(xj.)W(xj.)) 

A 

\ V ^ _ 2 



V| x ,Tr( 



Ax ± \ /_ Axj 

X -L 7T- )W [x±-\ — 



(52) 
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Fig. 5. The medium-induced gluon radiation diagram. 



with xj_ = (xj_ + x^_)/2. If the medium is large in the transverse direction, 
the dependence in x_|_ drops. 

The expression derived above is general. In the particular case of the 
multiple soft scattering, we find 

(pi) = -Ly dx + q(x + ) (53) 

Thus, we can interpret the q parameter as the momentum broadening per 
unit length [37] . 

Note also that within this approximation, in a homogeneous medium the 
number particle distribution at a certain transverse position after passing 
though a medium of lenght L is 

Af (x ± ) = J ^ e - ix ^^ (p_l) oc Tt(w*(0)W(xiS) (54) 

Thus, the particle distribution follows a diffusion equation in transverse 
space: 

d + M=^=qV 2 p± M. (55) 
From this fact, the interpretation of q as a jet transport parameter is clear. 

5. The medium induced gluon radiation 

With the formalism developed in the last sections we can now compute 
the medium-induced gluon radiation needed for jet quenching studies. In 
this formalism, the fast particle is produced at a given point inside the 
medium where the hard process, h, takes place, see Fig. [5j Then this 
particle and the emitted gluon suffer multiple scattering described by the 
path integral propagators ([26]) . We will work in the approximation p + S> 
k + S> k± but keep terms in k"j_/k+ as explained in the previous section. 
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Pi Pf 




(a) • (b) 

Fig. 6. Radiation diagrams after a hard process represented by the blob. 



5.1. The gluon emission vertex 

Let us first compute the amplitude for the radiation when the gluon 
does not interact after it is emitted. We start with the case when there is 
no scattering either of the quark, i.e. we consider the case of the radiation 
off a quark which has been produced in a hard process, h, with amplitude 
A4h ~ see Fig. [6]- we fix pi± = to have 

(56) 

where we have used the high-energy approximation to write 

p f -e pi-e k ± • e± 



Pf ■ k pi - k k_L 

The spectrum of radiated gluon is computed as the radiation over the 
elastic cross sections, where by elastic we mean the hard process A4h- In- 
troducing the relevant phase-space factors we find the double differential 
spectrum of radiated gluons in the vacuum 

all™ a s C R 1 
k+ dk + d^ ± ~ ^^kj (58j 

where Cr = Cf = {N 2 — l)/2N for quarks and Cr = Ca = N for gluons 
correspond to the average over the color states - see eq. Eq. ([55]) 

is the leading 1/z contribution to the vacuum radiation spectrum. For the 
case where one scattering is present - Fig. [6] (b) - the amplitude i^| 

~ -2g [ dx + dx ± igT b A b _{x)T a] ^± e - l{p f +k) ^u{p f ) (59) 
J kj_ 



6 For clarity, we omit the term Mu which cancels when the spectrum is computed - 
see eq. ((58} - if the momenta q are soft enough as we will standardly assume. 
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This expression contains a factorization of the radiation amplitude (|56p and 
the collision term which can be identified as the first term in the expansion 
of the Wilson line. The corresponding contribution resumming an arbitrary 
number of scatterings for a gluon emitted outside the medium is then 

^rad _ _ 2g J dx ± W{x±,x 0+ ,L + )T a k ^ ' f 1 - e- l ^f +k ^- x ^u{p f ) (60) 

Where xo+ and L + are the positions where the medium begins and ends in 
light-cone variables - so that the x+-integral in the Wilson line has these 
limits - and the subscript q is included to signal that the gluon does not 
interact with the medium. 

The case when the emitted gluon interacts with the medium can be com- 
puted similarly and the corresponding contribution from the vertex would 
give kij_/ki_i_ 2 , which is now an internal integration variable. In order to 
deal with it, we will write the delta function for the momentum conservation 
in the radiation vertex as 

(27r) 4 «5 4 (-^i " Pn+l +Pn) = I d A y S~ k -^ + ^-y. (61) 

and integrate over the three momenta k\, p n +i, Pn- By doing this, y signals 
now the spatial position of the radiation vertex. Let us define the contri- 
bution from a gluon radiated at y from a quark with momentum p n before 
the splitting and p n+ \ after it - see Fig. 

V g Kn) = J d*y e-^-^fl J ^VS)(*i) 

x exp {— i(kj + i — kj) ■ Xj} exp {— ik\ ■ y} . (62) 



-2 5 



k j+ A. 
k 2 - -He 



.( Xj ) 



7 We have used here the steps explained in Section 14.11 in particular eq. (|24[) , to 
simplify the triple gluon vertices in the high-energy limit 
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(we write fc n +i = kf for the final momentum of the gluon). Writing e_ = 
ki_i_ ■ej_/fti+ an d proceeding as explained in section [4~2l it is easy to see that 
the — component of (|62p . which will be used later, is 



9 



d y d ~x.\_\_dxi + e 



k + e± dy± 



G'o(y±,y+;xi_ L ,xi + ) 
igA-.(xi) Yl J d 2 x j ±dxj + G (xj- 1± ,x j ^ 1+ ;x j± ,x j+ )igA„(x j )e l]si -t' ± ' Xn± (63) 



3=2' 



so that the sum over scattering centers gives^l 

d 



V 9 ~ / d V e 



-i(p, 



n+l 



-Pn)-y 



dy± 



d x_ L G(y_ L ,y+;x n _ L ,x n+ ))e" 



And we define for future convenience 

;2 



i d 
dy ± 



ik-r i -x i 



d x 1 G(y i ,?/ + ;x ril ,x n+ ))e 4 ^ 



ik /± -x ± 



(64) 



(65) 



In the way eq. ([6211 has been defined, the amplitude for the process 
q — > qg corresponds to make the replacement A(y) — > A(y) for one of the 
scattering centers in the derivation of the quark Wilson line described in 
Section 14,11 so that the calculation is now similar to the one performed 
before. More specifically, the amplitude of Fig. can be written as 



M™ 



n 

i=0 



d A v 



t j4 



(2tt)< 



d Xi 



-tp-x 



■n 

3=0 



( 2gT""> ) ^ ^ Xj ' ' Pi (jPi-tei-vj+i) 



Pj + ie 



- n 

fc=n+2 



(-2 5 T afc ) 



^ ttfc (^fc)-Pfc cipfe .fa fe -^^ 



(frr<*"u(pf) (66) 



The structure in poles and phases is the same as the one studied in Section 
14.11 The poles give again an ordering in the longitudinal variable x+, so 
that the position of the radiation contribution A- {x + , xj_ ) is ordered with 
the rest of the external fields ^4_(x4-, xj_) to give 



M 



rad 



dx„ 



J dyxe i{pf ~ Pi)x ' yx Pexp ig J 
i2gT b A\x + ,y ± )Vexp 



Hi 



dy+A(y + ,y ± ) 



dy + A(y + ,y ± ) 



(67) 



We include here the case when the gluon does not interact once it is emitted insia 
the medium, otherwise a factor Go(y±, y+; x„i, x n +) should be subtracted. 
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+ 



+ 



Fig. 8. The three contributions to the squared amplitude of the medium-induced 
gluon radiation. The dashed line is the cut indicating the final outgoing particles. 



We have made explicit the color matrix T b at the radiation vertex while it 
is included as a redefinition of the external fields in the rest of the cases as 
done in previous sections. We are interested in the case where the quark is 
completely eikonal, so, we fix yj_ = to get 

M™ d = [ L+ dx + [ dxle- ik ^W{0;x 0+ ,x + ) x 

^+ Jx 0+ J 

xT b e ± ■ ^-G b (y x = 0, x+; x ± , L + )W(0; x + ,L + ) (68) 

The total amplitude for the medium-induced gluon radiation is then the 
sum of ([60D and §8$) 

M™ d = M™ d + . (69) 

We will now compute the spectrum of radiated gluons in the presence of 
a medium, including all the relevant color factors to perform the medium 
averages. 



5.2. The medium-induced gluon radiation 

The locality of the medium averages - see also below - allows for a simple 
diagramatical interpretation, in which three different cases appear when 
the amplitude (|69|) is squared depending on the position of the radiation 
vertex: when the gluon is emitted inside the medium in both amplitude and 
conjugate amplitude; when it is emitted inside the medium in amplitude and 
outside the medium in conjugate amplitude; and finally when the gluon is 
emitted outside the medium in both amplitude and conjugate amplitude - 
see Fig. [51 We take the case that x + < x + to obtain - see eqs. (|7Tj) - ([76"j) 



(\Ma->bc\ 2 ) 



9 

N 2 -l 



2Re 



1 

i? 



nr / dx j 



dxA 



dxdx e jk - 



_(x-x) 
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x+ x + 
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a 


ai 




b 

c 



Fig. 9. The color structure of the medium averages in eq. (|70| . 



/ c 9 

(W aai (0, x 0+ ,x + )T^ ibi —G ClC (y = 0,x + -,x,L + )W blb (0,x + ,L + ) 
2 k, f L + 



± [ + dx + [ d^ x (w aai (0,x 0+ ,x + )T^ ibi 



h , h 



X 



^G ClC (y = 0,x + ;x,L+) x W w (0, x+, L + )T b %wl 



+ 



+ ^ (70) 



Where we have written an eikonal Wilson line W(x±, xi + , X2+) for each 
quark propagation between positions x\ + — X2+, the gluon propagator^] by 
G(xj_, yj_, X2+) and the corresponding vertex factors as determined in 
the previous section. We have also explicitly included the color indices. 

All the medium averages can be written in terms of the Wilson loop 
average (|43p for gluons. To see how this works, it is useful to draw the 
conjugate amplitude as the Wilson line for the corresponding antiparticle 
in the amplitude - see Fig. [H 

Taking into account the composition of propagators with x + < z + < y + 

G(x±, x + ;y ± ,y + ) = J dz±G(x±,x + ;z_ L ,z + )G(z_ L ,z + ;y±,y + ) (71) 

One can make the substitutions 

W M1 (0,0,z + )w£ a (0,0,x + ) — > Wi iai (0,x+,x+) (72) 



9 Notice a small change of notation here on the order of the variables for the propagators 
G to more easily follow the lines in the figures. 
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W blb (0,x + ,L)wl- bi (0,x + ,L) -^W blb (0,x+,x + ) (73) 
Where we have used that in the averages 

'\V aai (0, 0, x + )wl(0, 0, x+)) = 5 ai - c (74) 
y c6 (0,5 + ,L)Wt i (0,5+,L)) = 5 c5l (75) 

The locality (in x+) of the medium averages we have presented in the pre- 
vious sections imply that the only non-zero contribution come from Wilson 
lines which overlap in x + . Eqs. (|72D and (I73D tell us that the contribution 
from the quark rescattering outside the longitudinal range cancel 
when doing the average. The remaining contribution from the quark lines 
can, in fact, be written as a Wilson line in the adjoint representation - i.e. 
as if the quark- antiquark propagation in the segment (x+,x+) would be a 
gluon propagation because - see the Appendix: 

TXwl- ai {Q,x + ,x + )Tll_W blbl {0,x + ,x + ) — > W c i- Cl (0,x + ,x + ) (76) 

With these simplifications, only possible in the case that the quark is taken 
as completely eikonal, eq. fjTOj) can be written in terms only of Wilson lines 
in the adjoint representation where the medium averages (|43p for gluons can 
be used. 



-±r / L+ dx+ [ L+ dx, I dxdxe*^*-*) 



xTr(w A (0, x+, x + ) — G(y = 0, x+; x, L + )— G(5c, L+; y = 0, x+] 



■ — -± y dx + y dxe ik ^ x Tr(^(0,x + ,L + )— G(y = 0,x + ;x,L+) 

, 4CW 2 



,(77) 



Using the relation eq. (I71D we have the average of only two Wilson lines at 
each segment in Fig. [9]- we write explicitly the two transverse components, 
then y and y will be put to again: 

Tr(iy A (y,x + ,x + )G'(y,x + ;x,x + )) (78) 
Tr (G(y, x + ; x, L+)G(x, L + ; y, x+)) (79) 
Tr(W A (y,x + ,L + )G(y,x+;x,L + )) (80) 
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The averages (ITSjl and (|79l) - taking x + = L + - give 



1 



N 2 -l 
T>r exp 



Tr (W A (y, x + ,x + )G(y, x + ;x, x H 



2 



exp 



1 



n(£)cr(r) 



(81) 



with the boundary conditions r(a;+) = y — y and r(x+) = x — y. The other 
average gives [42j HH [44] 

yTr (G(y, x + ; x, L+)G(x, L + ; y, = 
P+ (-2 



Vti I T>T2 exp < i / d£ 



y (*?(£) -*■!(£)) + 2<£Mri-r 2 ) 



Pu / Pv exp <^ i / 



^A(fl • v(£) + in(0a(v(0) 



(82) 



where we have made the change of variables 

u (£) = ri + r 2 

v(0 = ri - r 2 (83) 
The integration in u gives a (^-function constraining v to a straight line 



v s (e) = v (x + )^^ + v(* + )> ; 



— x + 



(84) 



the result is 142 



1 



N 2 -l 



Tr (G(y, x + ;x, L+)G(x, L+; y, x+)> 



(2^) 2eXP {S [(y " X)2_(5 '" i)2] "^ / WWO)} (85) 



with 



w(0 = (y-y)^— ^ + (x-x)^±; A£ = L + -z + (86) 



A£ A£ 
Eq. (fSTI) and ([86]) allow to compute the average of the first term in ((771) 

e ik x (x-x) 



dx<ix- 



Tr (VT A (y, x+, x+)G(y, x+; x, L+)G(y, x+; x, L+)) 



N 2 - 1 



j^—^f dxdxdze*^*-*') Tr(W A (y,x + ,x + )G(y,x + ;z,x + )) x 

xTr (G(z, x+; x, L + )G(x, L + ; y, x+)) = 
: /" dze ik ^- z /C(y-y, 3 ; + ;z + y-y,x + )exp -i / ^n(£>(z) (87) 
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Where we have defined 

K, (r(x + ), x + ; r(x+), x + ) = / Prexp 



LJx+ 



Putting all together, we obtain the radiation spectrum in the presence of a 
medium. 



k 



dl 



1 L L + «n(0<T(x) d d 



as ° R 2Re ( L+ dx + / rf 2 x r"*-* • 



dx+ e 



h + J x+ 



dy <9x 



K(y = 0,x + ;x,x + ) 



~ 2 TT " ^""^(y = °' x +! x ' L +) 



+ (89) 



5.3. The multiple soft scattering approximation 

In section 14.31 we have presented different average procedures to solve 
with (I88p . In the multiple soft scattering approximation, valid for 
opaque media, the dipole cross section is approximated by its quadratic 
term and the medium averages of two Wilson lines are given by (|43p . In this 
approximation, the path integrals <\88\) correspond to a harmonic oscillator 
of imaginary frequency 



/C (r(x_| 



T>r exp 



^1 di 



r 2 +i 



g(0 



(90) 

In order to proceed, we need to say something about the temporal depen- 
dence of the transport coefficient. Two classes of media have been studied: 
a static medium in which q(£) = q is a constant; an expanding medium in 
which the density of scattering centers is expected to produce a dilution 
as q(€) = 9o(?o/0 a > with a characterizing the speed of the dilution, and 
a = 1 for the Bjorken scaling scenario - see Section 12.11 Explicit solutions 
for this path integral and the corresponding spectrum (|89p are given in the 
Appendix. In the next sections we will present numerical calculations of 
these spectra, study their properties and explain how they are included in 
present phenomenology of jet quenching in heavy ion collisions. 



5.4- Numerical results and heuristic discussion: static medium 

In Fig. [10] we present the results for the double-differential medium- 
induced gluon radiation spectrum for a quark traversing a static medium. 
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Fig. 10. Left: numerical results for the medium induced gluon radiation spectrum 
of a quark in a static medium as a function of the dimensionless variables (p?Tj) . 
Right: Same but integrated in kt < lu 



The results are given as a function of the variables 

co c ^-qL 2 k 2 ^ (91) 
2 H qL y ' 

One important feature of the spectrum is the presence of small-/cj_ and 
large-u cuts which can be understood by the formation time of the gluon 



^form — i 9 • (92) 

kl 



The presence of this coherence length can be traced back to the non-eikonal 
terms in the propagators (|3ip . Recalling that these terms come from keeping 
k 2 _/2p + corrections in the phases and translating light-cone to ordinary 
variables 

^(..t-CWW) „ (93) 



n 



these contributions define the coherence time (|92p . When if orm <C L mul- 
tiple incoherent collisions are present and parametrically the spectrum is 
proportional to L/tf orm . In the opposite limit, when tf orm 3> L the gluon 
formation time is much larger than the medium size and the whole medium 
acts as a single scattering center. As a result, a reduction of the gluon ra- 
diation is produced in the last case. This is the generalization to QCD of 
the Landau- Pomeranchuk-Migdal effect |39|, BOl B5] . The numerical effect 
appears clearly in Fig. [10] as a suppression of the spectrum for small val- 
ues of k 2 . An important consequence is that the spectrum is not collinear 
divergent (i.e. it can be safely integrated to k± = 0) nor infrared divergent 
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(i.e. it can be integrated to w = 0) as can be seen in Fig. [JUJ In con- 
trast the vacuum part of the spectrum ([89]) present both collinear and soft 
divergences 

^ °* C * 1 1 (94) 



dujd?k\ it 2 k 2 



The position of the infrared cut in the /cj_-integrated spectrum - Fig. [TUl 
Right - can also be understood in terms of the formation time: integrating 
the spectrum in the kinematically allowed region < kt < m, and noticing 10 ! 
that (k 2 ) ~ \fquj then a suppression of the spectrum for to < q 1 ^ 3 should 
appear. 



5.5. Expanding medium 

In the physical situation of a heavy-ion collision, the medium formed 
is rapidly evolving and expanding both longitudinally and (as data seem 
to indicate) also transversely. A simple way of including at least part of 
the effect of this evolution in the formalism is to consider a medium which 
is diluting with time, so that the density of scattering centers decreases 
with some power-like behavior ra(£) ~ The corresponding transport 

coefficient will follow a similar behavior q(£) = <7o(£o/£)° as discussed before. 
In the Appendix the corresponding spectra for different values of a are 
computed. For practical applications, however, it is very convenient to use 
a scaling law which relates any expanding medium to an equivalent static 
scenario with time-averaged transport coefficient 

2 f L+ €° 

Q=tfj (95) 

The quality of this scaling for different values of a is presented in Fig. [TTJ 
Although some deviations appear for small values of R = uo c L they are, 
in practice, of limited relevance as they occur only when the radiation is 
anyway small - small lengths and/or q. 



5.6. Massive quarks 

It is now also simple to introduce the case of massive quarks into the 
formalism in the limit M <C E. In order to do so, let us observe that when 
solving the poles in pj_ in each quark propagator, eq. (|14p now reads 

r e i(x 1+ -x 2+ )k- 2vrz f M 2 1 

/ dk -^ — — = ~®{x2+ ~ xi+)- — exp^ (x 1+ -x 2 +)\. (96) 

J 2p + k^ + te 2p + { p + J 



This can be estimated by taking (fcj 2 ) ~ gfform, using eq. (|92[) (kt) ~ V<& 



lectures 'jq printed on February 2, 2008 



29 




Fig. 11. The gluon energy distribution for different values of the expansion param- 
eter a =0, 0.5, 1 and 1.5 and the corresponding scaling using Eq. ([55)1 . Fig. from 
[38]. 



This introduces only a multiplicative factor in the eikonal trajectories. This 
factor depends on the energy of the quark. The relevant +- component of 
the momentum is thus p+ before the gluon splitting and (1 — z)p+ after 
the gluon splitting. So, the total contribution from the amplitude times the 
complex conjugate amplitude to be included in the integrand of ([89]) is 

f M 2 \ ( M 2 

exp <^ % (x + - x+) > x exp <^ i- — (x+ - L + ) 

I P+ J I (1 - z)P+ 

( m 2 a r m 2 , 

exp < — i (xq+ — x+) > x exp < — i— — (x + — L + ) > (97) 



P+ J I (1 - Z )P+ 

which, expanding 1/(1 — z) — 1 — z and writing in terms of the gluon energy 
is just [36] 

f x 2 M 2 A . . 

exp<i— (x + -x + )} (98) 



So, in the limit of small-z and M <C E, the leading mass correction to the 
medium-induced gluon radiation corresponds to include the multiplicative 
factor flSgJ inside the integrand of ([85]) - see also [36l H71 B8lli9] . 

In Fig. [12] we compare the spectra of gluons radiated off massless and 
massive quarks. The effect of including quark masses into the propagators 
is to reduce the amount of radiation. This is known in the vacuum under 
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Fig. 12. fcj (left) and angular (right) spectra of medium-induced radiated gluons 
off a massless (red), charm (blue) and bottom (green) quarks [36] . 



the name of dead cone effect, as the mass terms in the propagators define 
a minimum angle Ode ~ M/E below which the radiation is suppressed. In 
the medium the situation is a bit more complicated because as we have 
discussed, formation time effects already suppress the radiation for small 
angles. In this case, two competing effects exist, on the one hand, there 
is a genuine suppression produced by mass terms in the propagators; on 
the other hand, the formation time is now smaller if or m ~ 2lo/(M 2 + kV) 
and the LPM suppression appears at smaller values of k±, enhancing the 
radiation at small angles as compared with the massless case. This effect 
is, however, restricted to a limited region of phase space, negligible for 
most observables, and the net effect of the mass terms is to reduce the 
amount of radiation. The energy loss of heavy quarks is, in this manner, 
smaller than that for light quarks and the formalism predicts an one-to- 
one correspondence between the two. This smaller radiation translates, in 
general, into smaller suppression of high-p_|_ heavy than light quarks. 



6. Application of the formalism for jet quenching studies 

As we discussed in Section O the cross section for producing a hadron 
h at high transverse momentum px and rapidity y can be written in the 
factorized form 



da 



AB^h 



dpi,dy 



dx2 dz 

2 2 Z 



xifi/A(xi, Q 2 )x 2 f j/B (x2, Q' 



*ki 



dt 



-^l mcd (^/4) 



(99) 
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where a summation over flavors is implicit - see e.g. [50] for an explicit 
expression including kinematical limits, etc. For proton-proton collisions, 
parton distribution functions, fi/ p (x,Q 2 ), and vacuum fragmentation func- 
tions, D™ h (z, fj/p), are known from global fits to experimental data using 
the DGLAP evolution equations [12] . In the nuclear case, the PDFs are 
known from similar global fits [22 [231 El ES ESI E3 EE]- On the other 
hand, the FF, D^f h (z, contain the information about the medium we 
want to study. In the following we explain how the formalism of medium- 
induced gluon radiation can be implemented in (|99p to compute the medium 
modification of the fragmentation functions and how to use them to char- 
acterize the medium properties. 



6.1. Medium- modified fragmentation functions 

The theoretical description of the fragmentation of a high-py particle 
in the presence of a medium is not completely known from first principle 
calculations and some degree of modeling is needed. One possibility is to 
include all the modifications of the fragmentation functions into a modi- 
fied splitting function in the DGLAP evolution equations. For a simplified 
discussion let us take into account only gluon FF - including other flavors 
translate into a summation over flavors - whose DGLAP evolution is 



d „ , . f 1 dz a s „ . . „ fx \ 



(100) 



where t = [i 2 F and P gg {z) is the splitting function describing the probability 
that a daughter gluon has been radiated from a parent gluon with fraction 
of momentum z. The probabilistic interpretation of the DGLAP evolution 
is more clear from its equivalent (at LO) integral formuation (see e.g. |13j ) 



D(x,t) 



A( (W s, W + A W f f ^ / £ P(z)D (« tl ) . (101) 



The first term on the right-hand side in this expression corresponds to the 
contribution with no splittings between to and t while the second one gives 
the evolution when some finite amount of radiation is present. The evolution 
is controlled by the Sudakov form factors 



A(i) = exp 



to 



(102) 



with the interpretation of the probability of no resolvable branching between 
the two scales t and to- The definition of the Sudakov form factors and 
its probabilistic interpretation depend on the cancellation of the different 
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divergencies appearing in the corresponding Feynman diagrams, see e.g. 
|51j . Although such cancellation has never been proved on general grounds 
for partons re-scattering in a medium, it has been found in [52] that, under 
certain assumptions, all the medium effects can be included in a redefinition 
of the splitting function 

P to \z) = P vac (z) + AP(z,t), (103) 



This possibility was exploited in |54j where the additional term in the split- 
ting probability is just taken from the medium-induced gluon radiation by 
comparing the leading contribution in the vacuum case - see eq. ([3]) 

, , 2irt dl med , . 

AP(z,t)~— — — , 104 
a s dzdt 

Implementing fllQ3f> and (|104p into (llOip the medium-modified fragmenta- 
tion functions can be computed - see also [55J for a related approach. Only 
the initial conditions of the evolution need to be specified. In [55] the KKP 
set of FF [56j was used for the vacuum as well as for the medium at the initial 
scale Ql i.e. D meA (x,Ql) = D vac (x,Q%). In this model all the medium- 
effects are built during the evolution. The motivation for this ansatz is the 
following: in hadronic collisions, particles produced at high enough trans- 
verse momentum hadronize outside the medium. So, this assumes that the 
non-perturbative hadronization is not modified by the medium, whose ef- 
fect is only to modify the perturbative associated radiation^!. All present 
radiative energy loss formalisms rely on this assumption. 

The fact that the medium-induced gluon radiation is infrared and collinear 
finite allows for a simplification of this formulation, valid when E ~ Q ^> 1. 
Under these conditions, eqs. (jlOip and (|102p can be written as 



D™ d (x,t) =p D™(x,t) + J J±-p(e)D™ (r^r e ,t) (105) 



with the probability of energy loss given by a Poisson distribution 

dl med 



Po = exp 



did) I dk±- 



(106) 



p(AE) = p ^ [ f[du>i ^ dk 

k=i ' J U=i ^° 



dl med (cui) 



dwdk 



o 7 dudkj 

sl^Uj-AE ](107) 



11 For processes with different kinematic conditions [57] this assumption could not hold, 
but the negligible effects seen in dAu data at RHIC [10] indicate that this is a rea- 
sonable assumption for particle production at high-pt in nuclear collisions. 
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R z 

Fig. 13. Left: The probability that a parton loses no energy when traveling through 

the medium (|106|) as a function of the variable R = oj c L — ^qL 3 - Fig. from [35]. 
Right: Medium-modified fragmentation functions computed in a DGLAP approach 
(54] and compared with the quenching weights results explained in the text. 



and e = AE/E. The total distribution is normally written as 

P(AE) = p 6(AE) +p{AE) (108) 

Eqs. (|106p . (|107p and (|108|) are normally called quenching weights. They 
have been first proposed in [53] and together with eq. (|105j) constitute the 
basis of most of the phenomenology of jet quenching in heavy-ion collisions. 
Notice that the probability that a parton loses no energy when traveling 
through the medium, given by pq, can be large if the medium length and/or 
the transport coefficient are small. In Fig. [T3]we plot this probability and 
the corresponding medium-modified fragmentation functions. Also plotted 
are the results from [54] when the medium effects are included at every 
individual splitting - notice that for the case of hadronic collisions, the 
most relevant virtuality is precisely that of the order of the initial parton 
energy, fj, F ~ E. 

A similar formalism is used for heavy quarks, for which, as we have 
seen, the energy loss is predicted to be smaller due to suppressed radiation 
controlled by mass terms (|98p . 

6.2. The suppression of the inclusive high-px yields 

Fig. [T3] shows clearly that one of the effects of the additional medium- 
induced radiation on particles produced at high transverse momentum is 
the softening of the spectrum due to the energy loss of the leading particle. 
The typical values of the fraction of momentum in hadronic collisions from 
(|99p are z ~ 0.5 — 0.7. The softening of the FF translates, in this way, in 
a suppression in the production of particles at high-pr in the presence of a 
medium. 
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Fig. 14. Left: Nuclear modification factor, Raa, for light hadrons in central AuAu 
collisions [58] , Data from [5U]. Right: Raa for non-photonic electrons with the 
corresponding uncertainty from the perturbative benchmark on the relative b/c 
contribution [61] . Data from [62] [63] 



In the massless case, the only parameters in the medium-induced gluon 
radiation are the transport coefficient and the length of the traversed medium 
The later is given by the geometry of the system while the formed is a fitting 
parameter - all the medium properties are encoded in q, so that measuring 
it we will learn about the properties of the medium. Although the geometry 
could, at first sight, seem a trivial feature in the calculation, different ge- 
ometries lead to different values of the extracted properties of the medium. 
In Fig. O the suppression computed with a static medium and medium 
density given by a Wood-Saxon parametrization is presented. The fit to the 
light meson suppression leads to quite large values of the transport coeffi- 
cient [581 EH] 

g~ 5 ... 15GeV 2 /mi (109) 

but determined with a large uncertainty. This feature can be understood 
as due to the dominance of surface emission: those particles produced close 
to the surface have also large probability to exit the medium essentially 
unaffected - see po in Fig. [13] In order to reproduce the large suppression 
observed, the value of q needs to be large, but at some point increasing this 
value translates only into a small reduction of the skin from which the par- 
ticles abandon the medium unaffected. It's worth mentioning here that the 
use of more sophisticated medium profiles, as given e.g. by hydrodynamical 
simulations, could lead to slightly smaller results [65J. 

In the massive case, once the value of q is known, the formalism pro- 
vides a prediction with no extra free parameter. At present, the suppression 
for charm or bottom mesons has not been directly measured. These quan- 
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tity is indirectly measured through the non-photonic electrons, expected to 
be mainly originated from semileptonic decays of the corresponding heavy 
mesons. In Fig. [II] the suppression from charm and beauty is presented. 
Also presented is the combination to the final electron yield from charm and 
beauty decays as given by the Fixed-Order-Next-to-Leading-Log (FONLL) 
description of the heavy quark production [63] • The band is just the effect 
of varying the masses m c and m\> as well as the corresponding renormaliza- 
tion scales as standardly uncertainties are computed in perturbative QCD 
calculations. The description of the data is reasonable but an experimen- 
tal separation between the charm and beauty contributions will lead to a 
definite answer on whether other effects could be present. RHIC experi- 
ments will be upgraded to address this issue. The LHC experiments have 
capabilities to distinguish charm and bottom - see some predictions in Fig. 
ESI 
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Fig. 15. Raa for D's (left) and for electrons coming from bottom decays (right) at 
y = for 10% PbPb collisions at y/s = 5.5 TeV/A, for different q (in GeV 2 /fm). 
Fig. from [66 . 



6.3. Jets in heavy-ion collisions 

The amount of information about the medium extracted by measuring 
inclusive particle suppression is limited as we have just seen. More dif- 
ferential observables could lead to a more accurate determination of the 
transport coefficient and to pindown the interplay with geometry. Relevant 
questions are where and how the energy of the initial parton is lost by in- 
teractions with the medium. For that, reconstructing the whole history of 
the parton while traveling through the medium would clearly be an ideal 
situation. Fortunately, some experimental handle on this is provided by jet 
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measurements. 

Jets are extended objects composed of collimated bunches of particles 
contained into a reduced region of phase space extending in azimuthal an- 
gle, A<1 >, and rapid ity, Ay. A jet can be characterized by the opening angle 
R = ^JA§ 2 + Ar] 2 defining a cone around the jet axis where the energy is 
deposited, but more sophisticated definitions are possible [68J. The descrip- 
tion of jets in the vacuum, as e.g. in experiments of e + e~ annihilation, is 
one of the most precise tests of perturbative QCD. With different degrees 
of refinement, a jet is computed as the output of a parton branching pro- 
cess starting by a highly virtual quark or gluon produced in the elementary 
partonic collision. Each branching is controlled by the splitting function, 
known at LO, NLO and NNLO. 

In the medium case, with modified splitting probabilities, the jet prop- 
erties are modified. We have already presented modifications to the longi- 
tudinal structure as given by the fragmentation functions. Jet properties 
are, however, richer and the transverse structure is expected to present a 
broadening produced by the multiple scattering in the medium. How the 
parton branching process is modified by the presence of a medium is still not 
known and, again, some modelization is needed. Clearly, the modification 
of the splitting function f)103|) provides a suitable model for implementa- 
tion in, e.g. Monte Carlo generators. Simple estimates of the effect were 
performed in Ref. [67] where the jet energy redistribution as given by the 
medium-induced gluon radiation (|89|) was studied. 

One of the main experimental issues in jet measurements is the calibra- 
tion of the jet energy - all the energy from particles belonging to the same 
jet need to be reconstructed. Two sources of error appear with opposite 
dependence on the definition of the jet size R: i) out-of-cone fluctuations, 
which decrease for larger jet sizes, due to radiation happening to angles 
larger than R ; ii) background fluctuations, which increase for larger R, 
consist on energy from other sources which enters the jet cone. Clearly, 
the high multiplicity environment of a heavy ion collision complicates the 
background subtraction which is under study by the different collaboration 
at the LHC [HHl EE EH E2] . Noticing that most of the particles in the back- 
ground are produced at small transverse momentum, the simplest way of 
controlling them is by introducing small-pT cuts in the measurement. 

At RHIC, real jet measurements in AuAu collisions are not possible due 
to the limited kinematic reach in transverse energy of the jet and the associ- 
ated large background. For this reason, two- and three-particle correlation 
analysis are performed as an alternative way of studying the jet structure. 
Here, one "trigger" particle is chosen within some range of p^ lg and corre- 
lated with "associated" particles at a relative angle A$ in a range of p|? soc 
usually smaller than the one for the trigger particle. A typical hard process, 
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with two back-to-back partons in the transverse plane, will appear in the 
final state as a two gaussian-like structures one centered at A$ = 0, corre- 
sponding to the trigger particle, and one at A<I> = tt, corresponding to the 
parton in the opposite direction. 

An important step forward is the first measurement of two particle az- 
imuthal correlations at large transverse momentum, with negligible combi- 
natorial background [73]. By choosing 8 < p^ lg < 15 GeV the correlations 
of associated particles with p^ soc > 6 GeV present the typical two-jet event 
topology both in pp, dAu and central AuAu collisions. These data support 
the picture of a very opaque medium with large energy losses, but with 
a broadening of the associated soft radiation hidden underneath the cut- 
off. Lowering this transverse momentum cut-off needs of a good control on 
the background subtraction, but the different collaborations agree in the 
presence of non-trivial angular structures [73]: the two-particle-correlation 
signal around the direction opposite to the trigger particle presents a dip in 
central collisions, in striking contrast with the typical Gaussian-like shape 
in proton-proton or peripheral AuAu collisions. The origin of this struc- 
ture is still unclear. In the presence of an ordering variable (as virtuality 
or angular ordering in the vacuum parton shower) the implementation of 
(|103p with the usual Sudakov form factors (|102p produce similar angular 
structures for energies uj < 2g 1//3 ~ 3 GeV in central AuAu [75]. At present, 
the most widely accepted interpretation is in terms of shock wave formation 
due to the hydrodynamical behavior of the energy deposition in the medium 
[76] . Cherenkov gluons have also been proposed as a possible source of these 
non-trivial angular structures [77J. Theoretical improvements on the used 
models, together with new experimental data, would help elucidating the 
origin of this interesting effect, which takes place, however, in a rather 
complicate kinematical region between the soft and the hard parts of the 
spectra. 

6.4- Is there an energy dependence on q? 

In all the application we have discussed q has been treated as a con- 
stant. This is a good approximation in the description of high transverse 
momentum particles at RHIC. However, the LHC will explore jets in a much 
broader range of energies which could allow to study the energy dependence 
of the momentum broadening. 

As we have seen in subsection l4.3[ q can be related to saturation physics. 
In fact, the momentum broadening of a probe propagating through the 
plasma can be expressed as 



(110) 
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Fig. 16. Jet quenching parameter q as a function of the jet energy for different 
path lengths. The square (triangle) marks the the value of q for thermal particle at 
T = 0.4 GeV (T = 0.6 GeV). L c = 1/xT is the coherence length, i. e. the length 
bellow which the particle scatters coherently with the medium. Above this length, 
no dependence on L is observed. Figure from |93) . 

where R is the color representation of the probe Q 2 S is the saturation scale of 
the medium and L is the length. The quenching spectrum depends on the 
gluon broadening. The saturation scale depends on the value of x probed 
in the interaction of the particle with the medium 

In a thermal medium, the saturation scale for typical thermal particles is 
given by /Ud ~ gT [93]. As the energy of the probe increases, the saturation 
scale increases because of evolution. In [93] it is argued that the typical x 
probed by a parton propagating in the medium is 



where E is the energy of the probe. Thus, the saturation scale grows and 
so does q. Since saturation is a non-linear process, it is not surprising that 
a non linear dependence of the Q 2 S with the medium length appears also. 

This effect has been studied in [93], where a simplified treatment of the 
evolution, based on the leading logarithmic approximation was performed. 
Their findings are summarized in Fig [161 A significant dependence on the 
jet energy of q is observed which could be measured in the LHC. As also 
proposed in [93], the broadening could be measured directly by studying 
the jet acoplanarity as a function of the jet energy, which constitutes a new 
observable for heavy ion collisions. 
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7. New developments 

In Section [5] we have seen that the spectrum of radiated gluon can be 
expressed in terms of the expectation value of two adjoint lightlike Wilson 
lines. The phenomenological applications of the formalism described in Sec- 
tion [6] are based in a particular model of this object, based on perturbation 
theory. As we have seen, the value of q extracted from the data is 3 — 4 times 
larger than what one expects from perturbation theory. This fact suggests 
that strong coupling non-perturbative effects are important. 

The computation of these effects in QCD is hard. Lattice QCD gives 
non-perturbative access to equilibrium properties. However, it is very com- 
plicated to extract dynamical quantities from the lattice, since the compu- 
tations are performed in euclidean time. What is more, the object we are 
interested in involves lightlike trajectories, which is an intrinsic Minkowski 
property. 

This difficulty has lead recently to significant interest in strong coupling 
computation in N = 4 Super-Symmetric Yang Mills in the infinity number 
of colors limit. This is, like QCD, a SU(N C ) gauge theory. However, unlike 
QCD, it is a supersymetic theory with no conformal anomaly. It also posses 
a different matter content, since apart from the common gluons, it contains 
fermions and scalar in the adjoint color representation. Despite of all these 
differences, computations in this theory are of interest due to a recent the- 
oretical development, the AdS/CFT correspondence. This correspondence 
allows to address the dynamical (real time) properties of the gauge theory at 
strong coupling in a regime in which traditional non perturbative methods, 
such as lattice, are not applicable. 

Even though the relation between QCD and M = 4 SYM is still question- 
able there are experimental (strong elliptic flow and particle suppression) 
and theoretical indications (lattice computations of the coupling constant) 
that suggest that the relevant physics of QCD close to the phase transition 
are those of strong coupling. Thus, the computations within the framework 
described below are of great value, since the AdS/CFT correspondence is the 
only method known to address these problems in a dynamical gauge theory. 
One of the great success of this approach is the computation of rj/s = 1/4tt. 
Current fits of hydrodynamical models to the elliptic flow yield a value of the 
shear viscosity close to that predicted from the AdS/CFT correspondence. 
This is a remarkable achievement for a theory in principle so different from 
QCD. Below we will review another computation in which a similar suc- 
cess has been obtained, the jet quenching parameter, and the momentum 
diffusion coefficient. 
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7.1. Brief introduction to the AdS/CFT correspondence 

In this subsection we will introduce the basic concepts of the AdS/CFT 
correspondence [78] necessary for the computation of lightlike Wislon lines 
within this framework. It is not the scope of this lectures to give a systematic 
introduction to this subject, which can be found in excellent reviews such 



In its weakest form, the AdS/CFT correspondence states that type II B 
Supergravity (SUGRA) in a particular gravitational background, the Anti 
de Sitter space times a five dimensional sphere AdS§ x S$, is dual to the 
large number of colors (N c — > oo), strong coupling limit of N = 4 SYM, 
a conformal field theory. SUGRA is a supersymmetric theory of gravity in 
10 dimension and it can be obtained as the low energy effective action of 
Type II B string theory. From this point of view, SUGRA is the effective 
theory that describes the string excitation with energy E 2 « 1/a' , and a' 
the string tension. 

The AdS$ space is a five dimensional space with constant (negative) 
curvature, R . A convenient choice for the AdS^ x S§ background metric 
is given by 



ds 2 = (-dt 2 + dx\ + dx\ + dxfj + — dr 2 + R 2 dn 5 , (112) 



where df^s is the metric of the five dimensional sphere. As can be seen, 
in the boundary, r — ► 00, the metric tends to the usual Minkowsky metric 
(times a prefactor). 

The duality maps the operators of the field theory to the boundary 
values of the SUGRA fields. In the gravity theory, all the dynamics are de- 
scribed by the classical equations of motion of those fields with the bound- 
ary conditions dictated by the correspondence. The expectation value of the 
operators in the field theory are obtained by evaluating the action of the 
gravity field on the classical solution. The coupling constant of the gauge 
theory, gyMi is related to the string tension and the curvature as 



where we have defined the 't Hooft coupling A = gy M N c . The SUGRA 
approximation is valid if R 2 » a' , thus, the 't Hooft coupling A — > 00. 

In order to describe the field theory at finite temperature, it is necessary 
to introduce a "black brane". From the SUGRA point of view, this object 
is analog to an extended black hole which leads to a 3 dimensional horizon 
at a given value of the coordinate r$. In the presence of this object, the 



as [79] [80]. 




R 2 




(113) 
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metric Eq. (11120 is modified 

ds 2 = (-fdt 2 + dx\ + dx\ + dxl) + %\dr 2 + R 2 dU 5 , (114) 

where / = 1 — r 4 /r 4 . The Hawking temperature of the black brane is the 
temperature of the field theory. In the metric Eq. (I114p . the temperature is 
given by 




7.2. Wilson Lines in AdS/CFT 

According to the correspondence, an external heavy probe in the fun- 
damental representation of the gauge group ("quark") is described in the 
dual theory by a classical string which pends from the boundary of the AdS 
space |81| . The string is a two dimensional surface parametrized by the 
internal coordinates r. a F^l 

(r,a)^X M = (t(T, C 7),x i (r, ( 7),r(T,a)) , (116) 

where the index i = 1,2,3. In analogy with the relativistic point particle, 
the action of the string is given by its area; this is known as the the Nambu- 
Goto action 

Sng = -^7 J dadr^ -det {G M N<9 a X M c> b X N } , (117) 

where a,b = r, a and Gmn is the background metric Eq. (| 1 12[) (Eq. (|114j) 
if we are interested in a finite temperature system) 

As we have obtained in the previous sections, the propagation of a heavy 
probe in the field theory is described in terms of a Wilson line in the fun- 
damental representation along the path of the quark. This path is also the 
boundary condition to solve the equations of motion derived from Eq. (|117p . 
After the solution of the classical equation of motion is found X¥ , the ex- 
pectation value of the corresponding Wilson line is given by 

(W) oce~ iSNG ( x ^ (118) 

As a simple example, let us consider a single massive quark at rest. Its 
space time trajectory is describe by a line at constant position xq. Due to 



The string also has coordinates in the 5-dimensional sphere. However in all the appli- 
cations below these remain constant and, thus, we will suppress them for simplicity. 
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symmetry, the solution of the Nambu-Goto action with the lowest action is 
the surface 

J# = M,r) • (119) 

According to the correspondence, the expectation value of the Wilson line 
is 

f i fT/ 2 ffmax ] 

(W) (x exp { -— / dt dr\ , (120) 



2na' 



-r/2 



where r max is an upper cut off. At zero temperature we can set tq = 0. In 
this case, the divergent prefactor is naturally interpreted as the mass of the 
quark 

2ira' v ' 

At finite temperature, the free energy of a single quark is given by 

Fi = — = (122 

2-na 1 2 v ; 



It has been also shown in [82] that as long as M » —F\ the thermal mass 
of the quark is given by 

M th = M — — ^— (123) 

7.3. Computing the Jet- Quenching parameter 

We now turn our attention to the computation of two parallel Wilson 
lines in the fundamental representation [831 

Tr 



W^(x 1 = -|)w c (x=|)]) , (124) 



where C is the line of constant rapidity 

tc = t cosh rj , X3C = —t sinh rj . (125) 

Lightlike trajectories are obtained by taking the limit n — > oo The line cover 
the interval t' = (— T'/2,T'/2) with T >> y. In this limit, we can assume 
that the dynamics are independent of t' . Due to this, it is convenient to 
introduce the following change of coordinate (boost to the rest frame) 

dt = dt' cosh rj — dx' s sinh n (126) 
dx^ = —dt' sinh n + dx' 3 cosh r\ (127) 
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To continue further, we need to choose a particular set of coordinates 
(r, a). Since the string solution we are interested on is independent of t', 
we choose r = t; it is also convenient to chose a = x\. It is also clear that 
the action is smaller if the the string world-sheet does not extend in the 
longitudinal direction x' 3 . Thus, we need to find a solution 

X% =(t', 3i,0,0,r(a;i)), (128) 

which should be symmetric under the change x\ — > —x\. Thus, if the string 
connects the two endpoints, it should have turning point at x\ = 0, i. e. 

r( Xl = 0) = 0, (129) 

where r = dr/dx\. 

The Nambu-Goto action reads, 



Note that this expression remains real only if r > tq cosh rj (r > tq since 
the string is outside the horizon). As we have seen, the maximum value of 
r , r max is proportional to the mass of the probe Eq. (|12ip . Thus, for a given 
value of the mass, there is a critical value of the rapidity 

M 2 , , 

coshr/ c = (131) 

W 

above which, the Wilson loop is exponentially suppressed. It is clear, then, 
that the limits M — > oo and rj — > oo do not commute. According to |84j . 
the lightlike Wilson loops are obtained by taking the rj — > oo first and then 
M — > oo. In this limit, the action is simplified to 

i 




S = 7 T'cosh?7 / dx i-^\m+— J ■ ( 132 ) 



Since the action does not depend on x±, we have a conserved Hamiltonian 

dC 

h = C-r- r , (133) 
or 

where 



E 2 r 2 
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From this expression, we obtain 



j) (1 - R 2 h 2 )) 



(135) 



Thus, in order to have real coordinates, R 2 h 2 << 1. Note also that in this 
approximation, the only point where the r vanishes in the horizon r = Tq. 
Thus, the string world-sheet always touches the horizon. 

From Eq. fll35|) we obtain an expression for the position of the string 



x\(r) 



R*h 



dr 



1 



Vl - R 2 h 2 

The constant of motion is determined by the condition 

x\{r — ► oo) 



y 

2 ' 



which leads to 



y 

2 



R 3 h 



r(l; 



r Vl - i? 2 /i 2 

Substituting Eq. (|136p into the action Eq. (|132|) we find 



(136) 



(137) 



(138) 



iV / ATv / vr^f--£ i 



r(I) \ 



i + 



r(|) 



7T- 



(139) 



where L = T' cosh 77 is the length travelled by the parton. 
At small transverse distance we obtain 



TV 




W c x 




(140) 



The term independent of the transverse position y corresponds to two in- 
dependent quarks. In [M] it is argued that this contribution should be 
subtracted. The remaining piece has the form of Eq. (|43p and we can read 



IF 



v / AT : 



!7 r 3/2 r(f; 



r(f) 



(141) 
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The index F indicates that the calculation has been performed with Wilson 
lines in the fundamental representation. This is not quite the object that 
appears in the radiative spectrum, where the Wilson lines were evaluated 
in the adjoint representation. However, these two objects are related in 
SU(N C ) theories [83] 

TrW A = TrWpTrWp - 1 (142) 
and, thus, in the large A^ c limit we obtain, 

q=^N c T\y^ (143) 

Despite of the fact that the above computation is performed in the in- 
finite coupling, infinite number of colours limit, we would like to have a 
numeric estimate as a comparison to perturbative expectations. To do so, 
we will assume typical values for a s = 0.3 and N c = 3 of RHIC physics. 
Recalling that a s = g 2 /Air, we obtain 

q = 4.5, 10.6, 20.7GeV 2 /fm for T = 300, 400, 500 MeV (144) 

Remarkably, the values obtained are much larger than the typical scale of 
the medium, set by the temperature. This is also what it has been observed 
in the phenomenological fits presented in Section [6j What is more, these 
numeric values are comparable to what was obtained in the phenomenolog- 
ical fits to RHIC data. Once again, this is a remarkable success, since the 
theory in which this parameter has been obtained, N = 4 SYM, is not QCD. 
This observation, together with the small values of the viscosity at RHIC, 
constitute a compelling indication (but, clearly, not an evidence) that the 
relevant physics at RHIC are those of strong coupling. 

However, we should be careful in comparing this values to RHIC phe- 
nomenology. As we can see, the jet quenching parameter depends explicitly 
in the number of colors and the coupling. This is different from other trans- 
port parameters computed in AdS/CFT such as the shear viscosity. The 
explicit dependence on the number of colors hints that, unlike the ij/s ratio, 
the result Eq. (|143p is not universal. This statement has been corroborated 
by explicit calculation in other theories with gravity duals. In fact, the 
number of colors is the only measurement of the medium density in the 
M = 4 SYM theory. Based on this and on the explicit computations in 
other theories, the authors of [M] have proposed the following scaling 



QQCD = \ qN=i ~ 0.6qaT=4 

V «Af=4 



(145) 
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where sqcd and sjv"=4 are the entropy density of QCD and M = 4 SYM 
respectively. Even with this correction, the obtained values are comparable 
with those obtained from the phenomenological fits. 



7.4- Momentum Broadening in Gauge Theories 
In the previous discussion, we have elevated the dipole formula Eq. (j43 



to a non perturbative definition of the jet quenching parameter, the medium 
scale responsible for the emission of hard gluons. The values obtained for 
this parameter are in surprisingly good agreement with those obtained from 
experimental data. This is a great success, however the correspondence is 
a powerful theoretical tool that allows to address issues that go beyond the 
estimation of the value of parameters. 

One of those issues is the relation of q with the momentum broadening. 
In subsection we have shown that in the dipole approximation q coincides 
with the mean momentum transferred squared per unit length. Since this 
is a quantity that could potentially be explicitly measured, it is interesting 
to check whether this relation still holds at strong coupling. 

The goal is to compute the momentum broadening of a probe within the 
AdS/CFT correspondence. We study a very heavy (almost) static probe. In 
this case, its motion in the plasma can be described by a classical equation 
of motion under a random force 

*-* < 146 > 

where, the force T is due to the thermal fluctuations of the (chromo) electric 
fields in the bath [85J 

T = J d 3 xQ\t,x)T a Q( y t,x)E a {t,x), (147) 

with Q(x,t) the heavy quark field. 

When the quark mass is very heavy, the time required to change its 
momentum is much greater than the typical scale of the medium, T. Thus, 
the distribution of "kicks" can be characterized by its second moment, 

{Fi(t)F0)) =K5 ij 5(t-t'). (148) 

where the constant k is the mean squared momentum transferred to the 
probe. From this, it is clear that 

dt(F(t)F(0)) HQ , (149) 
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where the average is taken with the partition function of the gauge theory 
plus the heavy quark. 

In this expression there is a technical subtlety. The correlator in Eq. (|149p 
is not time ordered and, thus, its extraction from the partition function 
should be treated with care. It is convenient to express k in terms of the 
retarded correlator. By using the Kubo-Martin-Schwinger (KMS) [86J rela- 
tions we can relate it to the retarded correlator 

2T 

k = lim ImG R (u) , (150) 

where, 

iG R (t) = e(t)([T(t),F(0)])HQ • (151) 

This is a rather technical step but it is motivated by the fact that there is a 
well defined prescription to compute retarded correlators in the AdS/CFT 
correspondence. 

The partition function of the bath in the presence of a heavy quark is 
given by [87] 



s 

= f d 3 xJ2(s'\Q(x,- T K PH QHx,-T)\s') , (152) 

s' 

are states of the gauge theory without the heavy quark. 
As is well know, the Hamiltonian in the partition function can be un- 
derstood as evolution in the imaginary time direction. The description of 
dynamical processes requires the introduction of a time contour that spans 
the real-time axis and extends towards the imaginary direction. This is 
known as the Schwinger-Keldysh contour which is illustrated in Fig. [TTJ 
In equilibrium, the partition function is independent of the initial time T. 
However, in computing a real time correlator it is necessary to specify the 
time in which the bath was created, since the phase of the heavy quark 
depends on it. 

As in the case of the energetic probe, the propagator of the heavy quarks 
field leads to a Wilson line. After some simple algebra, the froce-force 
correlator can be expressed as 

(T C [Htc)Hm H Q = ^- x 

(ix[U(-T - t c ) E(t c ) U(t c , 0) E(0) 17(0, -T)] ) , (153) 
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Fig. 17. Schwinger-Keldysh contour. The crosses represent the points of insertions 
of the electric field operators while the dashed lines are the paths of the Wilson 
lines connecting those insertions. The periodic boudary condition in imaginary 
time makes the force- force correlator Eq. (|153j) gauge invariant. 

where Tq denotes order along the time contour. 

The expression above can be obtained from a Wilson line along the time 
contour with a deformation in the y direction 5y 



Wc[5y] = T c exp 



dt (A + by A 



c 



y) ( ' 



(154) 



by taking (functional) derivatives on the path. 



{T C [T{t c )nm 



HQ 



S 2 W c [Sy] \ 

(W C [0]) \Sy(tc) Sy(0) / 



Sy=0 



(155) 



Thus, the broadening is obtained by studying small fluctuations of the Wil- 
son line along the heavy quark path. 



7.5. AdS/CFT Computation of the broadening 

As we saw in Subsection 17.21 a static Wilson line in the gauge theory 
is described in the gravity theory by a single string pending down from 
the boundary of AdS. Our task is to study how small fluctuations in the 
boundary along a path Sy propagate into the bulk. This is illustrated in 
Fig. [13 
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8y 



Boundary 



Boundary 



r=cosh T| r Q 



r=r o Horizon o Horizon 

Fig. 18. Left: Static string solution. The dashed line represents the fluctuations in 
the transverse coordinates Right: Trailing string correspoding to a probe moving 
at finite velocity. The string approaches logarithmically to the horizon. There is a 
new scale appearing at r — \J cosh r\ ro 



The string solution is a surface parametrized by 
Xf { = (t,r,0,y(t,r),0) , 



(156) 



which minimizes the Nambu-Gotto action, Eq. (|117j) . and with boundary 
condition 



y(t,r -> oo) = Sy(t) 



(157) 



Since the boundary conditions are set at r — > oo it is convenient to introduce 
the variable 



We also introduce the dimensionless variables 



TTTt 

irTy 



t 

y 



(158) 



(159) 
(160) 



Since we are interested in small fluctuations, we expand the action to 
quadratic order and we obtain 



Sng 



E 2 



dt du 



2-no 1 / 2u 3 / 2 



(161) 



where / = 1 — u 2 . Performing a Fourier transform 



y{t,u) 



. , . dtxi 

<■ ■— i/(»)y»(«) — i 



-imt 



(162) 
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where Y m (u = 0) = 1 and, thus, y(to) is the Fourier transformed of the path 
on the boundary. The equations of motion are 

^-^^^ + ^ = 0. (163) 

This differential equation has a regular singular point at the horizon, u = 1. 
Close to this point its solutions behave as 

y ro = (l-n) ±iro/4 F ± (n), (164) 

where the functions F±{u) are regular. These two solutions correspond to 
modes that are in- falling (-) and outgoing at the horizon (+). This can be 
seen from the fact that, close to the horizon, the lines of constant phase 
follow the path t = ±ln(l — u)/4. Retarded correlators are obtained by 
imposing the physically intuitive condition that all the modes close to the 
horizon are in- falling |881 [89] . 

The above prescription for the retarded correlators imposes the bound- 
ary condition at the horizon. Our task is to find the regular function F- (u). 
Inserting Eq. (|164p into Eq. (|163p we find 

2u — l- 3u 2 4 - u 9 \ „ . , . 

Since F_(u) is regular, we can expand it in a power series in to. Taking into 
account the boundary condition Eq. (|162j) we find 

F_(«) = 1 + 

— {^-4tan- 1 ( v / ^) -61og2 + 41og(l + y/u) + 21og(l+u)} + 

o 

0(tv 2 ). (166) 

Finally, the retarded correlator can be computed from the quadratic 
action. Following the prescription for Lorentzian correlators [88] we obtain 

G R (u) = -A(u)Y^(u)d u Y u {u)\ u ^ , A(u) = f . (1 67 ) 

Inserting Eq. (|166p into Eq. (|167p and Eq. (|150p we have 

— 2T 

K = lim IrnGflfw) , (168) 

= V\T 3 tt, (169) 
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where, as before, R 2 /a' = y/X. 

Our computation of the momentum diffusion constant k has the same 
temperature and coupling dependence than the jet-quenching parameter. 
This is not surprising, since the theory is conformal and the vA dependence 
is common in the observables that involve semiclassical strings in AdS. How- 
ever, the numeric coefficient is different. In fact, the relation between q and 
the momentum broadening derived in Eq. (|55p would imply 

q F = 2k , (170) 

which does not hold. However, before drawing any conclusion we must point 
out that the computation we just performed is in the static limit, while the 
q is in the light-like limit. Thus, in the next subsection, we extend the above 
computation to quarks moving with a finite velocity. 



7.6. Drag Force and Broadening at Finite Velocity 

We study now a single heavy quark propagating in the thermal bath at 
a finite rapidity rj along the x\ direction [82, 90J. As before, the quark is 
described in the gravity theory by a semiclassical string. Assuming that the 
quark is in a stationary situation, the string may be parametrized as 

X M = (t,vt + £(r),0,0,r) (171) 

with v = tanhr/ and £(r — > oo) = 0. With this parametrization, the Nambu- 
Gotto action is 




Since the action only depends on the solutions to the equations of mo- 
tions compatible with the specified boundary conditions are parametrized 
in terms of a constant of motion 

r 4 1 

— ./g 7 ; = C (173) 

R 4 llj cosh 2 7,-rS/r 4 , r 4 ,„ V ; 

V 7 + bF^ 

For these solutions, the classical action is given by 

Scl = —, r dr , f~ C2Ri / r ] (174) 
27TC*' J ro f (1/ cosh 2 rj - r 4 /r 4 ) ' V ' 

Note that the action has an imaginary part unless 

r 2 

C s = ±^smhn. (175) 
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Thus, for sufficiently long times, all other solutions are (exponentially) sup- 
pressed and the only relevant ones are those fulfilling Eq. (|175p . Inserting 
Eq. (fT75l) into Eq. ffT73|) we have 



E>2 2 

£' = ±^>, (176) 



which can be easily solved 



£ = ± - — ^arctan (—^ ~ arctanh (177) 

From these two solutions, only the (+) one is physical, since the string trails 
back form the quark EE This string solution is illustrated in Fig. [I8j 

As opposed to the static case discussed in the previous section, the 
sting bends behind the probe. Since the string has tension, work must be 
performed in order to keep the stationary situation. The momentum flux 
from the boundary to the bulk along the string is given by [82] 190] 



where X' , X represent the derivatives with respect to a and r respectively 
and 

-h=( K X'x)-(X') 2 (x)\ (179) 

the determinant of the induced metric. Substituting the solution Eq. (|177|) 
into this expression we find that there is a non vanishing flux of momentum 
in the x direction. This means that the probe loses momentum at a rate 
[821 [90] 



dp 1 TT 



TT 



-V\T 2 -== . (180) 



dt x 2 vT 

The reader may immediately find that the velocity dependence is that of 
the relativistic momentum. Multiplying and dividing by the mass we find 

f = -„ D p, v> = W)a»^ s . (i8i) 

At small momentum, the hard probe is expected to follow Langevine dynam- 
ics |91] where the average momentum loss has the drag form of Eq. (|18ip . 



A more rigorous derivation of this fact can be found in [82 , 90 
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For slow probes, the fluctuation dissipation theorem relates the value of the 
drag force and the momentum broadening 

ri D = K, — - — . (182) 
' 2TM v ; 

From this expression we see that the drag coefficient and the momentum 
broadening computed before fulfill this relation. The combination of these 
two independent computations can be taken as an explicit check of the 
AdS/CFT correspondence. 

Remarkably, Eq. (|18ip is, in principle, valid for all momentum. As in 
the previous section, we would like to compute the (transverse) momentum 
broadening of the probe at a finite velocity. As in the previous case, we 
need to study small transverse fluctuation of the string and compute the 
retarded correlator. The computation is greatly simplified if we introduce 
the following change of coordinates 



i = ■ ( t H — arctan (— ~\ arctanh (— ~\ (183) 

^sh^j V 2 \r J 2 \rJ y J 

— — -y/coshr/ arctan (^J cosh 7/ — ^ + - \J cosh rj arctanh ^a/cosIit/ ~j) > 

(184) 



\/ COsh 7] 



In these coordinates, the induced metric in the string h a t> = GMNd a X M dbX N 
is, 

h ti = -^f(f), (185) 
R 2 1 

hff = -T# 2 77^, ( 186 ) 

Hr = 0, (187) 

where /(f) = 1 — r^/r 4 . This is formally the same as the induced metric in 
the static string studied in section 17.51 but replacing the (t, r) coordinates 
with (t, f). 

Since the induced metric has the same functional form, the equation 
of motion for transverse fluctuations is also the same as in the static case. 
However, there is a subtlety: from the point of view of the moving string the 
horizon is lifted from its original location tq up to a higher scale ro^/coshrj. 
This is the scale where the local speed of light in the AdS space coincides 
with the speed of propagation of the probe. Thus, for excitations originated 
in the boundary, this scale plays the role of the black hole. 
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To obtain the momentum broadening we repeat the same steps as be- 
fore. The only difference is that the in-falling boundary condition should 
be imposed at the lifted horizon. Since, as we have remarked, the equations 
in these coordinates are the same, we do not repeat the computation and 
write only the final result |92j . 

k t = \J cosh r]y/XT 3 ; vr . (188) 

The obtained transverse momentum broadening grows with the rapidity 
of the particle propagating in the medium. This is an unexpected result, 
since the jet quenching parameter q yields an energy independent result. As 
mentioned, these quantities coincide in perturbation theory but they do not 
seem to do at strong coupling. However, before drawing this conclusion we 
should note that the computation presented in this section is not valid for 
arbitrarily large values of the rapidity. 

Since the quark loses momentum, the stationary situation is only reached 
when there is an external source. This can be thought of as a non vanishing 
electric field £ coupled to the quark living on the brane. The magnitude of 
the electric field grows with the velocity of the quark 

However, the brane does not support arbitrary values of the electric field, 
since it becomes unstable under pair creation. This occurs at a critical value 
of the electric field 



M 2 

e c = Tx . (190) 

As a result, there is an upper limit on the rapidity of the probe described 
by the solution Eq. (|177D 

( M \ 2 

cosh7 ? < — J . (191) 

The limiting value of the rapidity coincides (parametrically) with the one 
encountered in the computation of the jet quenching parameter, Eq. (jl31|) . 
Thus, the computations based on the trailing string cover a different range 
of rapidities than the computation of the jet quenching parameter. 

Thus, we cannot yet resolve the issue of whether the broadening and the 
q are related beyond perturbation theory. However, it is remarkable that we 
have obtained an energy dependent result for "small" values of the rapidity, 
while the ultra-relativistic limit is rapidity independent, as indicated by q. 
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Having a consistent dynamical picture that is valid in both regimes would 
be desirable and could lead to new observables able to point out whether 
the matter created at RHIC is indeed strongly coupled. 

The computation of the broadening is also interesting on its own. Slowly 
moving quarks follow Langevine dynamics and the heavy quark diffusion 
constant obtained from our computations 



is comparable with the value obtained from fits to the single electron sup- 
pression and i>2 (originated from charm and bottom quarks) within a Langevine 
model [91j . The energy dependence of the broadening constitutes a new 
observable which could, in principle, be measured by measuring the jet 
acoplanarity with fully reconstructed jets |93j. 

As a final remark, we would like to mention that the drag forced ob- 
tained for the heavy quarks is limited to the same constraint in the value 
of the rapidity. At very large rapidities (close to the light-cone) the mass 
of the probe should be irrelevant for its dynamics. Even though an explicit 
calculation has not been performed, it is clear that the mass dependence 
obtained in Eq. (|181 1) is not valid in the ultra-relativistic limit. 



The formalism described in Sections 2] to [6] and others [52U95] have been 
very successful for the description of high-p-p data from RHIC, specially 
for the inclusive yields. We have already commented on some limitations 
which would need to be addressed for the applications to the LHC, where 
the reach in the px spectrum will be enhanced by more than one order of 
magnitude and other observables as jets will be available for the first time 
in heavy-ion collision experiments. 

The main conclusion from the phenomenological studies of jet quenching 
is the characterization of the produced medium as extremely dense. We 
have quoted the values from [581 [59] of Q — 5. ..15 GeV 2 /fm. A perturbative 
calculation of the transport coefficient [M] underestimates this result by a 
factor of at least 5. This seems to indicate that the typical cross sections 
in the medium averages described in Section H] are much larger than the 
perturbative estimates based on a free gas of quarks and gluons. 

Several results from RHIC point, in fact, to this direction, as the small 
viscosity needed to reproduce the soft bulk data. Most of these phenomeno- 
logically accessible quantities are dynamical parameters for which the lattice 
calculations at finite temperature are not well suited. For these reasons, a 




(192) 



8. Last comments 
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flourishing activity on the relation of these findings with theoretically com- 
putable quantities in String Theory by the AdS/CFT correspondence has 
recently started. Although still highly speculative, this connection is open- 
ing new ways of facing the challenges on the study of collective properties 
at the fundamental level. 
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9. Appendix A 

Let W(x + , xj_) be the Wilson line at xj_ in the fundamental represention. 
Thus, it verifies 



-id+W(x+,x ± ) = A a _(x + ,x ± )T a W{x + ,x ± ) 
We now define the object 

W b A c (x + ,x ± ) = llY [w\x + ,x ± )T b W(x+,x ± )T c 
Taking the derivative with respect to x + we find 



(193) 
(194) 



- id + W A (x + ,xj_) = A°l-^Tt [^t( x+)X± ) 



rpb rpa 



Thus 

Using that 
we find 



- id+W% (s+,x_l) = Atif had Wf (x+,x ± ) 
(T<\) hd = -if abd 



W(x + ,x ± )T c (195) 

(196) 
(197) 



- id+W A (x+,x_l) = A_T A W A (x+,x ± ) 



(198) 



Thus, it is the Wilson line in the adjoint. The normalization guaranties that 
W at the origin is the identity. 
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10. Appendix B 

For consistency with previous approaches, used, in particular in Section 
El we present here the analytical results in ordinary variables - this means 
making the changes u ~ k + / y/2, x ~ x + j\[2 



dl 



u> 



dud 2 k ± ~ (2tt) 2 uj 



2Rc 



L+x 

dx / d 2 x e - ik ^' x 



1 f L+X ° I rl. 



UJ 



x ^~ ■ J~^(y = °> x ; x ' x ) ~ 2 tt ■ ^"^(y = °> x ; x > L ) 

ay ax ay 



a s C R 1 



+ 



7T 



2 k 2 , 



(199) 



Notice that the definition of the transport coefficient, q, differs also by a 
factor V2 [83], so the path integral is [59ll40ll4Tl [38] 



f C(r(x),x;r(x),x\.) = j 



T>r exp 



i— I d£[i ■ +i 



— i 

2u 



(200) 



which corresponds to a 2-dimensional harmonic oscillator with time-dependent 
imaginary frequency 



n|(go) = _ t m = _ { qo_ (to 



and mass u>. The solution of 
K(r%, x; T2, x\lo 



can be written in the form 



exp [iS c i(r 1 ,x;r 2 ,x)} . 



2iri D(x,x) 

Here, the classical action S c \ in (|202|) takes the form 



S c i(r 1 ,x;r 2 ,x) = — 



r d (0 • ^r d (0 



(201) 



(202) 



(203) 



where the classical path r c i(£) satisfies the homogeneous differential equa- 
tion 

' d 2 n 2 a (tjo 

_de 



r c i(0 = 0, 



with initial conditions 

r d (x) = ri and r d (x) = r 2 . 
The fluctuation determinant D (£,£') in (|202|) satisfies 



(204) 
(205) 



d 2 ^(&) 



.^2 



= o, 



(206) 
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with initial conditions 



= and ±D(Z,Z% = z, = l 



(207) 



In practice, -D(£,£') is found by combining the two independent (scalar) 
solutions f u / a of (|204"D . 



£>(e,O=^*(/i(0/2(O-/2(0/i(O) . 



(208) 



and fixing the norm by the initial condition (|20T|) . The solution of (|200p 
can be written in terms of £)(£,£') and two £- and ^'-dependent variables 

Cl, C2, 



/C (r(x),x; r(x),x\u) 



2ttD(x,x) 



x exp 



-Hi) 



2D(x,x) 



(drl + c 2 r2 2 - 2ri • r 2 ) 



(209) 



The particular form of the coefficient variables c\ and c 2 depend on a. 
Most of the phenomenologically relevant values are for a < 2 for which the 
two independent solutions of the homgeneous differential equation (1204|) are 

H3- 



(210) 
(211) 

(212) 
(213) 



/i(0 = y/tlufatoafa)^ 



where /„ and -K^ are modified Bessel functions with argument 

1 



2-Q 

In terms of the variable 

z(0 = 2i/n«(&)C^ . 

the solution (|209p is given by [we use z = z(£), z' = z(£')) 

D(CO = (2 ^ o)) 2, (**T [lu{z)K v {^) - K v {z)I v {z')] , (214) 



ci = z (^-\ [l v - X (z)K v (J) + K u -i(z)I v (z')] , 
c 2 = z (-V [K v {z)I v -x^)+I v {z)K v . x ^)] . 



(215) 
(216) 
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We now write the three contributions from (11990 as 



dl 



dujd 2 ~ki it 2 



a *C F {h + h + h)=u C " 



med 



dujd 2 \i 



dP 



dujd 2 \i i 



(217) 



with, I§ = \/k\. I A and I5 define the medium-induced part J med which, 
using d09|) - (I2T61) leads to 



4uj 2 



2Re 



dx 



-4^4 



x exp 

J 5 = — 2Re 

2uj 



4 (D 4 - i A 4 Ba) 



L+x 



X() 



dx — ^ exp 
5? 



(£> 4 - ^ 4 L? 4 ) 2 (L> 4 - iA A B A f _ 

(218) 

k 2 1 

(219) 



4A 5 5 5 



where 

A A 
M 



— , J B 4 = ci(x,x), L> 4 = J 
2D(x,x) 2 

S5 = ci(L + x ,x) . 



L+xo 



2D(L + x ,x] 



d£n(£)a(r) (220) 
(221) 



In the case of a static medium, a = 0, the functions I±i/ 2 (z) and K± 1 / 2 (z) 
entering (I209j) have explicit expressions in terms of exponentials and the 
results are simply 

M = 2sin(Q(Lx)) ' Bi = COS{n{2 ~ X)) ' Di = \ n ° C{L ~ S) • 
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